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Surface Displacements in an Elastic Half-Space”) 
By H. A. Lane 


Die vorliegende Arbeit ist ein Bericht über die Resultate einer Untersuchung von elastischen ,„‚Grund- 
wellen“, die auftreten, wenn eine konzentrierte Last (mit einer Heavisideschen Zeitabhängigkeit) einen 
Druck auf die Oberfläche eines elastischen Halbraumes ausübt. 

Die Resultate bestehen aus 


1. Integralausdrücken für die einzeln genommenen Beiträge der P, 8, M und der Rayleighschen Grund- 
wellensysteme zu den Verschiebungen und den zugehörigen Potentialen; 


2. bestimmten Ausdrücken und Kurven für die Verschiebungen, die jede Wellenart erzeugt. Für den 
Spezialfallv = 1/4 reduzieren sich die gesamten Flächenverschiebungen auf die bekannten Formeln, die von 
Pekeris erhalten wurden. 


This is a report of results obtained in an investigation of elastic ground waves initiated by a point load 
(having a Heaviside time variation) acting on the surface of an elastic halfspace. T'he results consist of 


I. Integral expressions for displacements and associated potentials everywhere in the half-space for the 
separate contributions from the P, S, M and Raleigh surface wave systems. 


2. Explicit expressions and curves for the surface displacements for each type of wave. The total surface 
displacements (evaluated for v = 1/4) reduce to known expressions obtained by Pekeris. 


Hacroaman pa60Ta ABIAeTCA MORAANOM O Pe3yJIbTATaX, HOJIYYCHHBIX IPH NCcJIeNOBAHHH 
9AHACTMYHBIX ‚,OCHOBHBIX BOJIH‘‘, BO3HURAIIMX, KOTMA COCPeNOTOyeHHaAH HarTpysaka (c Xe- 
BNUCA ÜNOBOÜ 3ABHCUHMOCTBIO OT BPeMEHH) OKA3bIBACT NAaBJIeHME HA HOBEPXHOCTB HEKOTOPOTO, B 
HalleM CcJIyyae 3NACTH4YHOTO, IIOJIYIPOCTPAHCTBa. 

IlosıyyeHHINbIle pe3yJIBTaTbI COCTOAT U3: 

1. UHTeTpaubHbIX BbIP&;KeHMÜ Alf BBATbIX B OTAEJIBHOCTH NOJIeH, BHOCHMBIX BeJImymHaMmu P, 
S, M u cucTeMmaMmu OCHOBHbIX BOAH PeüJen B CMEINEHUA U COOTBETCBYIOIIME IIOTEHIMAJIEI; 

2. oNpenesnleHHbIX BbIPa?KeHNÜ U KPUBbIX AI CMeileHuf CO3NaBaeMbIX KAsKIbIM BUNOM BOJIH. 
Ina vacTHoro cayyaa v» = 1/4 Bce CMeIleHnsT HOBEPXHOCTeÜ CBONATCH RK U3BECTHBEIM POPMyIaM, 
HONy4YeHHbIX Ilekepucom. 


Nomenelature 
The following nomenclature is used in the paper: 

t — time %,ß — ineident and reflection angles of P, 
1, 2,9 — eylindrical coordinates N NE 
R,6,9 — spherical coordinates na Or sin Ir 2v 
u, Ww — displacements (cylindrical) ee 07  sma 2(1 N 
U,W — displacements (spherical) N(«, ß) — tan cos?2ß + tan Psin?2 
A, u — Lame’s constants 0 — Sr 
v = Po1sson’s ratio Ri 2 
4A, w — dilatation, rotation ” 
8,W — scalar, vector potentials n ” Crt 

PER? 2 3 | # 
C;, =\/ = dilatational (P) wave speec h — sin ß 

a M(v) —-1—80%2+8(3 — 2sin?y) v* 

je en ; — (16 cos? y) v® 

= EEE 2 rot 1(8 ve speed 05° Y 
= 4 rn + gi(i = 1to 3) = roots of M(v) 
| CR = surface wave speed @, — [a — cos (« — 0)] [a — cos (@ + 0)] 
| o,T —= normal, shear stresses es — [b au u vn D cos (ß-+ 6)] 

0) = frequency H, = cos 2 P sin? P/N %, | 
k = separation constant Nö, — sin2 ßsin P/N(«, ß) 


Introduetion 


Wave propagation in an extended homogeneous, isotropie elastice medium has been dis- 
cussed by Lam [1], Stores [2], Love [3] and many other investigators. It is well known that, 
when the medium is of bounded extent, the treatment of reflections at one or more free surfaces, 


*) Paper No. P-1498, The Rand Corporation, Santa Monica, California. 
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or the satisfaction of boundary conditions resulting from the application of external load consi- 
derably increases the difficulty of obtaining appropriate solutions of the equations of motion. 
Among the general class of three-dimensional wave-propagation problems in a bounded medium, 
that of the propagation of axially symmetrie wave fields is of considerable importance. 

An examination of both two-dimensional and axisymmetrie wave-propagation problems 
in the technical literature indicates a number of observations which underlie the choice of method 
of the present paper: 

1. Transform methods frequently lead to integral expressions for displacements and stresses 
which are difficult to evaluate explieitly-particularly near the center of disturbance where the 
methods of stationary phase or steepest descent have been required to furnish useful results. The 
discussions in the recent work of GAarvIn [9], JEFFREYS [6], Larwoop [7], Hure and Core [11], 
[12], and Mexcner [10] amply illustrate this observation. 

2. Some of the particular methods applied to two-dimensional problems are capable, in 
principle, of extension to axisymmetric-wave fields. Yet the relatively few specific examples of 
axisymmetric-wave problems in the literature reflect the fact that the extension, in practice, is 
difficult to carry out. 

These observations led to the conjecture that the method used by SAUTER possessed certain 
distinctive advantages which, if extended from two-dimensional wave fields to axisymmetric- 
wave fields, would eliminate the aforementioned difficulties. The method, which is described 
briefly below, is related to methods used by SUMMERFELD and KoTTLer [4] in studying electro- 
magnetic wave-diffraction problems. 

The surface displacements of the problem considered in this paper have been determined in 
closed from by PEkERrıs [21]. Prkerıs [24], [25]!) has also explicitly determined the surface dis- 
placements for the closely related problem of the buried seismie pulse. Other investigators of 
the buried pulse problem are CAGnIARD [8], Pınwey [22] and Dix [23]. The solution for the 
surface pulse is obtained by letting the depth of the source approach zero. 


The Method of Sauter 


SAUTER’s method consists offorming an aggregate of solutions ofthe wave equation — doubly 
infinite in frequency and separation constant. Transform methods are used to satisfy the 
boundary load conditions. The resulting expressions for the displacements can be integrated over 
the real domain of frequencies by choosing a suitable path in the complex plane of the separate 
constant. The real part of the separation constant satisfies the condition of reflection (SNELL’s 
law) while the imaginary part adjusts the phases of the individual wave in the manner required by 
the boundary conditions. This effects a separation of the displacements (now single integrals) 
into separate expressions which, because of singularities in the integrand, can be identified with 
domains in the medium having different characteristics. The contribution of each domain and the 

wave behavior at critical points 
Lyt near the surface and on the wave 

front can be studied. 
This method was applied by 
pr  NAUTER to solve the two-dimensio- 
nal problem of a fluid halfspace 
initially at rest subjected to an 
arbitrary surface disturbance q(z, ?) 
[14]. The extension to an elastie 
half-space under the two-dimen- 
sional surface loading o, = o(z, D), 
Tr, = T(x, f) constitutes SAUTER’sS 
second problem [15]. Search of the 
literature failed to disclose the ex- 
tension to axisymmetricfieldswhich 
SAUTER planned. The present paper 
is based on the same heuristic physi- 
cal arguments advanced by SAUTER. 


X X 


Fig. 1. Elastie. Half-Space Showing Wave Regions (Left Half of Fig 
and Nomenclature (Right Half of Figure) 


ure) 


Statement of the Problem 
a 2 5 { = v a verlical downward unit force is applied at the origin of cylindrical coordina- 
esr, zin the elastic half-space shown in Fig. 1. The time variation of the force is taken to be the 


1 : ; 
) Reference [25] contains a rather com 


RR en plete discussion of the effect of S, P and $ P (diffracted) waves 
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Heaviside unit function. The waves that are propagated from this epicenter consist of a P or 
dilatational wave traveling with a speed C, and an S or rotational wave traveling with a slower 
speed Cz. Such waves have hemispherical fronts at any time £. In addition, a Mac wave having 
the straight line shock front and tangent to the S wave envelope exists. The surface (RayLzicn) 
wave traveling at speed C; (as indicated in Fig. 1) will be discussed later. 


It will be shown that the separation of wave displacements leads to separate wave systems 
denoted by I, II, III. The displacements (I) existin the P,M, S regions behind the P front of 
Fig.1. The displacements (II) exist in the S region behind the hemispherical front of Fig.1. 
The displacements (III) exist only in the M region. The time of arrival of a wave of specified type 
at any point Q is characterized by the (dimensionless) arrival times 


D a=.CztR ai 
(ID b='Crt/R b>1.0 
(I) b= Crt/R 1>2b> cos (# —y) 


Mathematical Analysis 
The displacement vector for a point in an elastic medium can be written as 
u=grad®Ö+ cu W 
where the scalar potential © satisfies the equation 


0:D 108 &0 1 98 


2D — —— == 
De ee RR 


and the vector potential  satisfies the equation 


CP) 20) 1CM) 18%) e 
0z2 or? Tr Or CH ot? Paar) $ . 


Mr) = 


The scalar potential ® determines the dilatational or P wave of speed C, while the vector 
potential, which has only one component, determines the rotational (S) wave of speed C7. 


Formal solutions, harmonic in time, are 


TEL 
i ® — kz— ot 


D=de 7 Jukr), 


which represent standing cylindrical waves of frequency & and separation constant k. 
The solution is generalized by writing: 


ze Lg Va 


—o 0 


Ju(kr) dk do, 
..@)% 


> i( Sr Zn) 
> g(k,o)e 7 er 


The stresses at the surface (z = 0) are 


1G,:) = u [ [Gkyo) Jılkr) ei! dkdo, 


— 00 


8 


ort) = u [ | Fk 0) Jukr) ot dkdo, 
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where 
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N 2 Ren es 
ko) — 2ki Vz ach (& ine re) Ms 
| = 


2 9i 2 >> 
= © EL 2kil/- — RRglk, 0). 
F(k, o) (2 K) fik, ©) + ri] 2 g(k, ©) 
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Using both a complex FOURIER transform and a HANker, transform?), | 


a Hrn 
>) [e>) k Zn j IR 
F(k, ©) e [ j Imu o(r', t) dot dr Jo(k r’) r’dr’, s Wr j | on | 
where o(r’, f‘) is the normal load on the surface. j ver 


For the Heaviside unit loading 


© . € oo > Mh 
0) Ik) eier dr’ dt = lim | ZZ Iukr) ar | eier ar — —E hf 
[Jen san: dr a = (kr See 
— © 0 v : {N N Er 
SI —ik 
since the convergence of the improper time integral can be justified®). Now F(k, ©) = re 
and Eqs. (4) reduce to re 


oo 0 


= — Jı(kr) iR? emiet dk do ale 
u AmtuoNkk, o) ’ 


ers | 3 

Se | re $ 

IRRE - i> 3 

o-|| Er dk do | ] Re 

—m@0 . 3 5 

[ ] denoting the bracketed terms of Egs. (4). fr 
This aggregate of cylindrical waves is equivalent to a family of plane waves propagated in 2 

all direetions, each with an amplitude dependent on the direction?). + 
u 


The integrals of Eqs. (5) are complex. To obtain real displacements, 
real and the separation constant k complex. Since a function and its c 
real parts, the frequency limits 
between the limits, 0, ©. 


we take the frequency & 


onjugate have the same 
are modified by taking twice the real parts of the displacements 


?) Reference [17] pages 8 and 12, 
®) Ref. [17], p. 9 or Ref. [18], P.:8317. 
4) See Ref. [18], p. 361. 
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To complete the frequency integration, an a i i ion i 
\ ‘ - ppropriate path of integration in the complex k 
plane is replaced by equivalent paths C, in complex « and £ planes defined by Peg 
D°. IR 
Kir nie: Ge een o0or ß plane 
implying the law of refraction, 
sinysina=sinß, where siny=(C7z/C, . - . (6). 


To obtain waves that decay with depth the paths €, are selected 
so that the imaginary parts of cos a, cos ß are positives as required by 
the exponential members of Egs. (4)°). If Eq. (6) were limited to real 
values of «, ß, the symmetry of the problem suggests that 0 <a < 90°, 
0=P= %°. The extension to complex angles requires 0 < Re (a) < 90° 
and 0 = Re (f) < 90°. These requirements are satisfied by the paths 
C, shown in Fig. 2. 


Eliminating the separation constant k reduces N(k, o) to „, er El 
(o/Cr)* N(a, ß) cot a, where® Fr IR 
N(a, ) = [tan a cos®?2$ + tan Psin?2ß]. in the. Complex wo 
The displacements become”) 
- R e io nn 50 C05s 0 — 
22? 4 Cru=—Re I» | Be sin « sin 0) sin ß H,(a) e a | 
c,o n 
. 'R z 5 io I eoß cos d— 
+ Re eos sin 9) cosßsin PH,(ß)e (er Aa ) maß, 
c, 0 - i 
* 1 
f ; 2 sa cos d— 
2m!uCrw=—Re | cota dl sin «sin 9) sin B.H,(e) e er) do da 
c,o0 % 
£ R io ER sß cos d— 
— Re | Es BE sin ß sin 0) sin? 8 H,(ß) e (er r-) do dß 
‚67 J 
co 


The interpretation of Egqs. (7) is as follows: The x members are referred to the path C, in 
the a plane. Furthermore, using Egq. (6), the integrand of the x members is a function of a only. 
This is emphasized by writing H,(«) despite the explicit appearance of $ in the definition of H,. 
Furthermore we shall later use alternative paths rather than C, in the a plane. Again Eq. (6) may 
be used to transform these alternative paths into equivalent paths in a f plane quite indepen- 
dently of the # members of Egqs. (7). Similarly, the integrand of the $# members of Eaqs. (7) is, in 
view of Eq. (6), a function of $ only and is here referred to a path C, in the complex plane. 


The singularities®) of the integrands are such that the path C, may be deformed into the 
path C, in either x or ß plane. The path C, approaches the angle a = Bin the a plane or = din 
the £ plane where 9 is restricted to the open interval 0 << 90°. Displacements on the axis 
and on the surface (9 = ©, 90°, respectively) are discussed later. 


Along paths such as C, or C, in either the a or plane, the frequency integral convergences?) 


since 


1 (72 a nn _—— ry 


J eo Jn(s ©) do = = = (n — 0, 1) - 
: J 


\r? + & 


5) All square roots in the paper are positive (except where noted). 

6) The same expression arises in SAUTER’s two-dimensional problem. 

?) It is convenient to use spherical rather than eylindrical coordinates for the point @ of Fig. 1. 
8) See later sections. 

9) Ref. [19], p. 64. 
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ir i | Oin the a plane and 
Th rgence requirement Re (r + is) 2 0 reduces to Im cos («+ > Oin 
im ce 0> On the plane. Both conditions are satisfied by paths C, und G,. The dis- 


placements reduce to 


—_ 08 (cos ß cos 8 — b) H,(ß) cos  dß 
223 Rusin 0 = — Im | A) da — Im ne 
6, j GC (8), 
sin? $ H,(P) d 
Bar Im [ara A) am [A 6 rt ] 


Ca 


G, = [a— cos (a — 6)] [a— cos (a + 9)], 
G, = [b — cos (ß — 0)] [b — cos (P + 9)]- 
The additional integrals which result from the frequency integration of the horizontal 
displacements u are discarded because they are real. They have the form 


Im [ H,(a) d« and Im J H,(ß) cos ß dß. 
[07 "2 


To see that these are real, it is sufficient to evaluate them across a path C, at infinity formed by 

‚closing the contour in the & and ß planes of Fig. 2. We may show that the integrals (from 0 — i oo 

to (”/2)— i 00) Im [ H,(a) da and Im [ H,(ß) cos # dß vanish. Since there are no singularities 
© Ö 


within the a or ß planes, by Caucay’s integral theorem the same integrals vanish on C,. 

The steps leading to Eqs. (8) depend upon the fact that the branch points k = + w/C, 
and k= + w/Cr of Egq. (4) are transferred ta= + /2, « = sin! (1/siny) in the « plane, 
respectively, and toß = + y, f = + z/2in the f plane, respectively, (Fig. 2). The branch cuts 
may be drawn from the branch points to + ooin such a manner that they lie outside of the regions 
of interest in the a, ö planes of Fig. 2. For this reason they are excluded from Fig. 2 and sub- 
sequent figures. The singularities (simple poles) corresponding to the zeros of N(a, ß) are discussed 
later. 


Separation of Displacements 


The separation of displacements depends on the branch points of G,, G, which result from 
the frequency integration. When a< 1.0, the quadratic factor in the « members of Eggs. (8), 
produces the branch points P, and P, (the branch cut is the heavy line shown in Fig. 3) correspon- 


a or 3'plane 


a orß plane 


Fig. 3. The en at 0, — Cy (u s 1.0 in the « plane) Fig. 4. The Closed Contour O3, 3,45 (for «>1.0inthe a plane) 
b =1.0,0<yin the ß Plane) (for 5 21.0, 0 < y in the $ plane) 


ding to « ah = cos"! a, For O <9d< 90°, applying Caucay’s integral theorem to the closed 
contour of Fig. 3, and discarding the real contribution from the path (, at infinity leads to the 
En the a Be of Eqs. (8) vanish when a < 1.0. This result corresponds to satisfying 
‚he radıation condition since outgoing waves can have no influence rior to reachi i 
} R aves ce B achin 
distance Q (Fig. 1) from the epicenter. h De 
When a > 1.0, the branch points are 0: Q, (ti i i i 
;’ 0, 8 ‚oınts are Q,, Q, (the branch cut is the heavy line) shown in 
Fig. 4) corresponding toa = + icosh-!a. The contribution around Q, vanishes and the result 
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When b < 1.0, the branch points are P, and P, (Fig. 3) and the ß members vanish by the 
argument used above. When 5b > 1.0 and 0 > y, the branch points are Q,, Q, corresponding to 
ß = 0 + cosh”! b and the ß members of Eqs. (8) again reduce to the contour CH 


- Corresponding to the phenomena of total reflection, the above arguments must be modified 
when 0 > y since diffracted P waves now arise from incident S waves. Physical reasoning suggests 
that along a path in the complex ß plane such as the horizontal path C, of Fig. 2, where Pis (very 
nearly) real, the exclusion of portions of the path corresponding’ 
to y< Re(ß) < 0 would, in virtue of Egq. (6), exclude imaginary 
contributions arising from the complex angle «. This suggests 
that the contour C, must be modified to C, as shown in Bier5. 20 
Then adopting the same arguments as before, it is inferred that 

-the appropriate limits are y to 0 — cos-! b for the M region and 
yto6—i cosh”!b for the S region when b > 1.0 and > y. An 
alternative argument is as follows: The branch cut can be extended 
from the branch point $ =y to Q when b > 1.0 and 9 >y. The 
final contour C, is the straight line Iying along this branch (shown 
dotted in Fig. 5). The limits and path for the M region is established 
by reducing b from values > 1.0 to values < 1.0. In this limiting 
prodedure Q, ascends vertically, coincides with 9 when b = 1.0 and ° 
is relabeled P, when b< 1.0. This -equivalent choice of path C, 
does not effect theresults (Eqs. (9)) which follow. a 


The displacements of Eqs. (8) may now be separated into the sets 
0 —icosh’ia 

1 "  a—cosacos9 H,(a) 

eng sin 0 Gi® 


Pplane 


N 
‚Alternative wQ, 
branch cut 


da, 


a0: 


8 —i cosh’ta 


cos 


0 —i cosh”!b 
—bH, 
u(II) = — Im ß ie «(P) dß | 


2? Ru sind G42 
). b>10, (9). 
6 —icosh’!b H (B) 
1 b 
w(lI) = — ER Im sin? ß er dp, | 
min (0, y) 
23 B cos 0 —b Hulp) 
1 cos ß cos d — 5 
u(1ll) = Sy Im | cos ß in GR dß, 
RR 1>b>cos(d—y), 
1 Berk) 
w(llI) = — ERn Im | sin? et dp 


The wave system I refers to the region behind the P-wave front of Fig. 1. The wave systems 
II, III apply to the S, M regions, respectively, in Fig. 1. For clarity, the paths of integration to 
be associated with Egs. (9) are redrawn in Figs. 6 and 7. 


Pplone 


a or plane 


9 2 
——— 


| 
ı 
| 
+ 
) 
I 
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Fig.7. Path for Waves (II) in 3 Plane, 6>y,5b21.0, 


Fig. 6. Path for Waves (I) in« Plane; Path for Waves (II) Path for Waves (ID) in $ Plane, 1 2b’ cos (d— 9) 


in ß Plane, 0 < y 
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Preliminary Discussion of Egs. (9) er R 
em (I) represents the (dilatational) P-wave. The rotation, w is proportional 


The wave syst 
to 52 (RU) +a z (RW) which vanishes. The scalar potential ® may be constructed from 
da 


the equalions 


od 1 08 
Ul-p Winde 


and is: 
0—icosh”!a 


1 VG de | 
ea j Al) 10a ein de: 


ß - 
This satisfies the wave equation, a result easily confirmed.by rewriting Eq. (1) in the form 


a 
lin) = 


Similarly the displacements II and III satisfy the requirement of vanishing dilatation 


(7) 10.70 : 
en ee —ı2 
(RU) 5 (RU) + nd 16 (W Rsin 6) 
The vector potential Y, obtained from the equations 
1 0) 3 1 0 S 
Un, == 29 Fr] (7 sın 0) Wir, sm = — Ban OR (R V sin 0) 
is 
„0Ö—i cosh’! ’H (B) 
1 A er 
use —_—, >10, 
u ano VGzdß, b> 
Min (6, y) 
0 — cos"!b ( 
Be Hp) \G >} > e= 
a ann VG, aß, 1>b>cos(d —y). 
2 
These satisfy the wave equation (2) rewritten in the more convenient form 
0? (7) 0° 
DENN : Be en : Se : a 

(b 1) Sp (7 sin 0) — cot 9 56 (7 sin 6) + 568 (Fsind)=0. 


The displacements and potentials vanish on the wave front. The boundary conditions on 

the surface 

5I+M)=0, HI D=0 for 1210; 
and 

5(I+1D)=0, e; 1I+1ID)=0 fr 12b>siny 
are also satisfied. 

The displacements (I) vanish on the wave front a = 1.0 in view of the limits of the integrals 
of Eqs. (9). Similarly, the displacements (II) vanish on the wave front b= 1 when d< y. The 
displacements (II) vanish on the wave front b = cos (d — y). The displacements II, III, must 
satisfy continuity requirements at b = 1,0 > y, but need not vanish individually. In this sense, 
b=1,9>yisnota wave front. These statements must be accepted provisionally at this stage!). 


The Surface Displacements 


lo examine Eqs. (9) at the surface, we use the substitutions 


D=sinp, ai. 2 Vsin? y — v (1 —20)—40v Y1 — v2 Ysin?y — v2 
N(a, ß) v M ; 
where M is the cubicin Y, (and Y — v2) 
« R 3 
1—-8Y-+8(@8—2sin2 v) Y?— 16 cos?y Y?] = — 16 cos? yI1l(w— 9) 
en ke 


10) A more rigorous proof will be given in sub i i 
u gorou: 8 subsequent papers. The proof depends upon integrating Eas. (9 
and passing to the limit as we approach each wave front. A similar procedure can be Be Oben ae 


ments on the axis. The arguments of this paper are not applicable at 8 = 0 because the branch points coalesce 
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having roots 9}, 93, 93. For any value of Poıssox’s ratio 0 <v 


that gj and g; are both less than sin? y while the root 9% > 2.0: 
Transforming the displacements (I) to an equivalent path in the ß plane (defined by Eg. (6)), 


< 1/2, the roots may be labeled so 


and using the substitutions 


v Vo? — sin? y(1— 2 vi)? 


lv) = 


MV\®—v | 
4 08 Vv® — 1 (v®? — sin? y) 
DvD) = Fe BD) 
92(0) MY» — » 
2 __ cın2 _— p2 u) 2 
ee sin?2y 1 — v2(1 En 


M Vb2 — v2 


> the surface displacements become 


ul, = Fr 2 v? r(v) dv 
sin y 
b 
b 


u(Ill) = Imeru (1 —2 v?) r(v) dv 


siny 
b 


1 : 
u) [E20 sinpy=b=T0, 
siny 
b 
w(III) = — — 2 v? g,(v) dv siny= Db= 1.08% 
n?ru 
0 (10). 
> 
HA ala 2 v? r(v) dv = 10, 
"10 
b 
— Dahl, 
u(ll) = Inru e 2 v2) r(v) dv = 
= » 
Br 2 — 202 d b>1.0 
0) [2 [ana Hz, 
sin y 1.0 
b 
— 1 
ex 2 er 2 d b=>1.0 
m(Il) = nr, 2 v2 g,(v) dv rg, P» g.(v) dv | 
siny 1.0 
The total displacements are 
b 
u(I+ en r(v) dv rzb 23%, 
== 
b 07 & 
er | r(v) dv b>10, 
= alle 
eg 
w(I+ 1) = g(v) dv r2b2Siny, 
3 ; 
w(I + ID) = er | gı(v) eg | 00 dv b> 1.0 


sin’ BD 1.0, 


sin we 1b 1.0; 


Fig. 8. Vertical Displacements I, II, I at Surface of Half-Space for » = 0.25 


Fig. 9. Horizontal Displacements T, IT, IIT at Surface of Half-Space for v = 0.25 
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When » = 1/4, siny= 1/3, M=8—4Y +56 ?— 32 y3 
identically to the surface en ren 5 ser 

The explieit results for Poıssox’s ratio, » = 1/4, obtained from Egs. (10) are given in the 
Appendix. The substitutions and integrals required to reduce Eqs. (10) to the results in the 
See are obtained from the procedure used by PEkErıs to determine total surface displace- 
ments. 

The surface displacements for » = 1/4 for each wave system, I, II, III (including the surface 
wave contribution discussed later) are plotted from the equations of the Appendix and are shown 
in Figs 8 and 9. The total surface displacements (denoted in the figures by I+ IlandI-+ III) 
were determined by Prkkrıs. The determination of individual rather than total surface wave 
displacements represents one advantage of the method adopted here. Other advantages are 
discussed in the concluding remarks of this paper. 


The Surface (Rayleigh) Wave 


There are no singularities on the path for the expressions in Eqs. (10) and (11) when b < 3. 
When b > g,, the pole at g, contributes additional terms. We may indent the path so it is semi- 
eircle about the pole 95 in the complex Y plane (see ImY 
Fig. 10). There is no real contribution to w(I) or w(I]). 5 
The contribution to u(I) and u(II) (for any Poısson’s 
ratio) is 


et er : 
16rrucoey (E—M(E— HM) Vr—g' 
b 1— 2.9) 1 
u(lI) = 32 2 2 n 2 . 2) VE 5 
arucos®y 9(5 — 9) (9 — N) Ver — g 
For » = 1/4, these reduce to REN re 
as + 5) 
LIE er 
16rru VB — g 
bIVY3 +1 
Ms Ei RN 
16rruV®— g 
—l 
re Ber 


sacru Vb: —_g 


CGoneluding Remarks 


The limiting surface displacements (of the Appendix) represent one example of the inter- 
pretation of Egs. (9). A second interpretation of an elementary nature can be obtained by con- 
sidering an acoustic medium (u — 0). Ifthe displacements everywhere (which arise only from the 
P waves) are differentiated with respect to time, the solution, corresponding to a pulse loading 
in an acoustic medium, is: 

l cos d \ {sin 0 
ruhe, an R- 


In front of and behind the pulse, the medium is unstressed and at rest. The pulse travels without 
dilatation, rotation, or acceleration as noted earlier by STOREs") and Love!?). Comparison of the 
two- and three-dimensional solutions exhibits the differences expected between wave propagation 
in even and odd dimensional space. In particular, only the two-dimensional pulse exhibits a 
residue or tail®). 

A third interpretation of Egs. (9) depends on explieitly completing the integration in order 
to determine whether the resultant displacements are capable of representing the complete wave 
field at all points of the halfspace. The discussion and contents of this paper constitute the first 
part of an investigation which was motivated by this third interpretation. As such, the paper 
limits itself to the integral relations and partial conclusions necessary to provide the background 
for the second part of the investigation. The second stage describes the procedure used to obtain 
the explieit interior displacement equations. There does not appear to be any reasonably coneise 
method of telescoping the entire investigation into a single paper. 

1) Ref. [2], p- 279. 


12) Ref. [3], p- 305. 
13) Ref. [20], p. 525 & Ref. [18], p. 388. 
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Appendix 


Explieit expressions for the surface displacements I, II, III are given below for Poıssox’s 
ratio = 0.25. The roots of M are 


Be: DES: Mr: 
N 1°’ 93 t 3 95 Te 


%. Pe 
32rurw(l): | 6 b2 (®?—H) 2+%(b: 2) 2 kl — b2) 2 


’ 


where: 


BeisBes. li. j  _VaVaı r 3 fa 

) : (11.18.— 19), k,=V3 3 +3 Vu + 19). 
32 a RL Be erh, 
32rurw(lll)=7-+ 652 5) "— gg) "+ kb! —)-2—k, (2 —b2) 2, 


wlm 


where: 


3 ; z 
k, = 'E lau 101, k, IE (11Y3 +19). 
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BO 

| ER 3 
32zurw)=—125 -6—2k(Y—b) ?, 32urw(ll) = 18 +12 —2k,(—b2) ?. 
Forb >9, 


Dr nu 6, 32rurw(ül) =18 4128. 


For u pl, 
Y3 


run) -— | yotiK -8E oma 4 (8-31) IIE20— 12 13) 1] 
45 +8 +33) I[(20 + 12 V3) R}} , 

Iararuti) = | BE 5K -omBR 43 WB) TIL — 1218) el 
+ (3 + 93) TI[@0 + 12 y3) re]}, 

R=,0®—). 
Korle=)b 0,5 a 

ar urul) -—/z oktıea— 309 (125?—1)K 97118] +33 3) [20-1273]. 
r +9 318) 1120 1233] 203). 
16a uruı) = — [3 kt E@R + 12m K a ms] + 8-13) 0 121%] 
+(3+ 3) nf20 +12 Y3]} . (16), 


2 
== 32 —1 
För.D —>.g,, 
u(l) = Eg. (12) + Egq. (15), u(II) = Eg. (13) + Eg. (16). 

The horizontal displacements contain complete elliptic integrals of all three kinds. The 
modulus k is specified below the displacements. The elliptic integrals are K, E, Il(n) where n is 
the parameter. The definition of the elliptic integral of the third kind is: 

n]2 
Il(n) = = > =” N er 
(1 + nsin? 9) /L,— K? sin? o, 
0 


Manuskripteingang: 24. 6. 1960 


Anschrift: H. A. Lang, The Rand Corporation, Aero-Astronauties Dept., 1700 Main St., Santa Monica, Cali- 
fornia, USA 
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Longitudinal Vortices in Shear Flow”) 
By G. S. S. LUDFORD 


Der Einfluß eines Längswirbels auf die Scherströmung an einer Wand wird auf Grund exakter Glei- 
chungen untersucht. Es zeigt sich, daß sich der Wirbeleinen Weg unter der Scherströmung und quer zu ihr 
bahnt und dabei diese von der Wand um einen Betrag abhebt, der von der Entfernung des Wirbels von der 
Wand, jedoch nicht von seiner Stärke abhängt. Es werden einige Vermutungen vorgebracht hinsichtlich des 
Verhaltens von Paaren und Gruppen von Längswirbeln. 


The effect of a trailing or longitudinal vortex on shear flow near a wall is investigated on the basis of exact 
equations. It is found that the vortex tunnels under and across the shear flow, Lifting the flow away from the 
wall by an amount depending on the distance of the vortex from the wall but not on its strengih. Some conjec- 
tures are made concerning the behavior of pairs and groups of longitudinal vortices. 


FIBJIeHne BONOUAMEerTOCA HAM IIPOAOABHOTO BUXPA B IIPOXOAAIIEeM BOJM3U CTEHKN TeYEHHN 
Nayuaercan Ha OCHOBe TOUHBIX ypaBHenni. OÖHapyskmBaetch, 4TO BUXpb Bblönpaer ceöe 
IIyTb ION NOTOKOM WIM IONeP&K Hero, OTPbIBAnH TeueHue OT CTEHKM B HEKOTOPOH Mepe, 
3aBucemmeli OT PacCcToaHNs BUHXPA OT CTEHKU HO He OT ero cmiIbI. Jlarıee TestamTch HeROTOPbIE 
IIPe1no.10:KeHns 0 TOBeNeHUN HAp U TPYIIH HPONONBHEIX BUXpeh. 


3) Research sponsored by the Office of Ordnance Research, U. S. Army, under Contract DA-36-034- 
ORD-1486. 
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$ 1. In a recent paper SCHULTZ-GRUNOW [1] has described various experimental conditions 


i itudi i i laminar or turbulent 
der which longitudinal vortices can be produced in the shear flow of a 
Hutare layer. The paper contains a bibliography, to which may be added the later work of 
Weske [2], with its appendix by BURGERS. See also GÖRTLER’s presentation [3]. 


Under suitable conditions (small enough speeds and smooth enough main flow) the motion 
is steady. Downstream from a disturbance element, such as a hemispherical cap, thin metal trip 
or a suction hole, there is a pair of vortex tubes which originate at or close to the element. For a 
considerable distance these tubes appear to follow the main flow direction, remaining remarkably 
parallel and distinct. In the case of a suction hole it is not diffieult to defleet away one of these 
tubes at its origin and so obtain, in an extensive region of the flow, a single concentrated vortex 


lying along the same general direction as the main flow. 


The object ofthe present paper is to give amathematical description of this last phenomenon 
and to draw some plausible conclusiong concerning pairs and groups of longitudinal vortices. For 
this purpose we shall consider an incompressible fluid and negleet the effects of viscosity, even 
though this is the agency which sets up the main shear flow. Observations indicate that its diffusive 


action is small. 


The author is indebted to Professor BURGERS for pointing out this problem. The basic 
equations were derived by BuURGERS in the appendix mentioned above. There is no difficulty in 
the case of uniform flow in the absence of a wall: a rec- 
2 U12) tilinear vortex may always be placed along the flow 
direction to form a combined motion. However, this 
cannot be done in the case of shear flow. We shall see 
that the presence of the wall is essential. The vortex is 
set at some angle d to the main flow (fig. 1) and it then 
has a tendency to move. This can be cancelled by the 
induced effect of the wall or, what is the same thing, of 
the image vortex, and thereby the angle ® is fixed?). 


Fig. 1. Rectilinear vortex in shear flow near wall Of course, an inclination 6 is also necessary for 
uniform flow in the presence ofa wall. However, certain 
features of the combined motion which are immediately evident for shear flow are obscured 
for uniform flow (due to its uniformity, for which displacement of the flow is not apparent). 
The vortex tunnels underneath the shear flow, lifting it away from the wall by an amount 
depending on its distance from the wall, but not on its strength. Hence the disturbance caused 
by an arbitrarily weak vortex is not of the same small order. Moreover, inside the tunnel 
the flow is to some extent arbitrary, presumably being fixed by conditions upstream at its point 
of formation. Thus as the wall recedes to infinity the original flow (uniform or shear) is pushed 
away to infinity, and whereas uniform flow along the direction of the now isolated vortex can be 
regained by suitable prescription within the tunnel at each stage, shear flow cannot. 


$ 2. We assume that in the original flow, undisturbed by the vortex, the velocity depends 
only on distance z from the wall and denote it by U(z). We further assume that the vortex stays 
at a constant distance z = a from the wall and at a constant angle 9 to the undisturbed flow, and 
ask for a disturbed flow which is independent of x, distance measured along the vortex. 


Steady flow patterns which do not vary in the x-direction satisfy the equations 


09V 9w ou au 
oy £ y a; 02 2 | 
N : 
dv ow 1 op ow ow 1 op ( 
D 2 | rn A dy® D- dy + w dr = — 3 Fe 


From these it follows that v, w are the velocity components of a plane flow in the y, z-plane and 
that u is constant on the streamlines of this flow: 


o 
u PR. 


Fr ne = UF. W).- un ee er 


D 


!) This is familiar also from lifting line theory, where each trailing vortex is actually inclined at a small 
downward angle to the uniform ineident flow, due to the downwash induced by the rest of the vortex system. 
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where y(y, z) is the corresponding streamfunction and fis as yet undetermined. The equation 
satisfied by y is, as usual, ? | 
022 


where the vorticity &(y) (in the plane flow) must be prescribed. Note that the vortex vector ofthe 
flow out of the plane is 


02 02 - 


ou ou 

I 5-5) 17@ 0») a ae 

so that the vortex lines are the streamlines of the plane flow, the vorticity on any line being 
proportional to the spacing between the lines: 


In the present problem the flow in the y, z-plane tends to the undisturbed shear flow 
v = U(z) sin 0, w = 0 at large distances, for which 


N ee 
0 


Hence, in that part of the flow which is covered by streamlines originating at infinity (y > 0, see 
fig. 2), o(y) is obtained by eliminating z between these two equations. Elsewhere it may be pres- 
cribed arbitrarily. The problem then is to find a solution of (3) for which 


v5 log Desert 20% 


=05 on! =, ee ED 


| 


vosind/UOd : as YPR2>o | 
0 


Here x, the strength of the vortex, is assumed 
to be given; the angle 0 must finally be chosen 
so that the vortex is free, i. e. 


V20 


Urz) sind 
ee 


0 7 
7 (#3, en)=0 ätı 7.0410): 


$3. This program can be carried out 
explicitly in the case of a uniform shear flow 


Ur z, 
for which 
o(Y) sin constant (7) Fig. 2. Flow in eross-planes for U=kz 
at points on streamlines which originate at infinity. For simplieity & will be given this value 

everywhere. Clearly the function 


1 Bar 2 
y= 5(k sin 0) z? + log n3,108 T; 


-satisfies the differential equation (3), when & is given by (7), and the conditions (5). The first 
term is the undisturbed shear flow in the plane; the second and third are the potential flows of the 
vortex and its image in z = 02). To satisfy (6) we must take 


R x 
RCRERER! 


. (8). 


Clearly there is a limit 
KEANUAR 


to the strength of the vortex for which this solution is valid. The streamfunction may now be 
written 
Za8 r 
(le 


2) In plane motion a uniform shear flow may be superimposed on any potential flow. 
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The pressure is given by e are ä 
P.< const. En +) —oy, 
0 2 
as is easily verified from (1). R 


Finally, there is the velocity component u out of the y, z-plane, which at large distances must 
tend to its undisturbed value (k cos 0) z. Since 


ee bi: 
Lan 7 % 


at large distances, we must take 


8a x 2\2 en 
Im = keoso]/ E ver) GORE a SE 


in (2) . This result only holds where the expression under the radical sign is positive, 1. e. at points 
on streamlines (p > 0) which originate at infinity. On the closed streamlines (y < 0) in fig. 2, 


u may be prescribed arbitrarily. 


On the streamline y = 0 which lies over these closed streamlines, u = 0 and one possibility 
is that u vanishes on them also. In this case the vortex produces a purely circulatory motion in a 
sort of tunnel underneath the raised shear flow. 


Equation (9) shows that all flows in the y, z-plane are similar. On the z-axis 
z\2 2% |2—a]| 
ae) HoelsFa| 


a graph of 2 y/x versus z/a is given in fig. 3. It shows that at its highest point the shear flow is 
raised to about twice (2.06) the height of the vortex above the plane. Note that this result is 


ee 
see 


’ 


ee l | 
vr 0 1 ? 3 4 5 
Eu u/ak c050 


Fig. 3. Variation of y on zaxis for U = kz Fig. 4. Variation of u on z-axis for U=kz 


ndependent of the strength of the vortex: the disturbance caused by a very weak vortex is not 
of the same order of smallness. Complete flows are similar only if the parameter x/a? k (i. e. 0) 
is.the same, see (10). 


As was mentioned before [see equation (4)], the vortex lines of the shear flow out of the y, z- 
plane are the streamlines of the flow in this plane, the vorticity on any line being proportional to 
the spacing between the lines. Thus the vortieity is everywhere increased by the lifting, and in 
RN Ne on the lowest line 9 = 0 [since f’(0) = 00, see (10)]. This is illustrated in 

Ig. 4, which gives the variation of u/a k cos O on the z-axis. The straisht line i is fi - 
presents the undisturbed flow. e SE ea 


BE This infinite veloeity gradient is not typical; U must vanish at the wall for it lo occur 
Consider the uniform shear flow 


U- k(z + 20) 3 
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which, with z, > 0, would perhaps be a better model for boundary layers. Equations (8), (9), and 


(10) are replaced by 


EN 


(9), 


sind = s 
4rna(a+2z,)k 
_.x*[?+2z2 rn, 
v ah a(a+2,) 1708 5 
Ip) = k.cos oz mw, 


eG 


= 2nVa(a-+ z,) 


% t 1) (z + 20)” 


Duehn, 
a ion. 


Now f’(0) is finite; also the flows in the y, z-plane have a similarity parameter & = a/z,, and the 
tunnel is closed at the points y= + yA x +3 (fig. 5). 


$4. The flow due to a single isolated vortex may always be superimposed on a uniform 
flow by placing the axis of the vortex along the direction of the stream. The fact that this was no 


longer possible for a uniform shear flow 
was the starting point of the present in- 


vestigation. The reason for this is now . 


clear. If in these last three equations we 
let @— co, keeping a+2z, and z+2% 
fixed, then certainly 6 —0. Butin the 
y, z-plane the flow tends to that of the 
given vortex, i. e. the shear flow is split 
apart by the vortex and pushed off to 
infinity in all directions. Whether this 
disruptive effect of longitudinal vortices 
on shear flow has a physical significance 


iz 


kzZ+Z,)siaO 


DELL, GEF EL 
oV4@+3 J 


Fig. 5. Flow in cross-planes for U=k(2 + 2,), 
2, positive (x = a/z,) 


is an interesting open question. 

In the solution presented above certain choices were made for the values of w and u on the 
closed streamlines. In fact they will be determined by the mechanism which produces the vortex. 
In WESKE’s experiments using a suction hole there is evidence that u is non-zero. This does not 
alter the above analysis, but only fig. 4. A different distribution of » would affect the solution, 
though probably only qualitatively. However, it leads to a difficult mathematical problem of the 
free-boundary type for the differential equation (3). In this connection there would be a case for 
eliminating the concentration of vorticity at the vortex and distributing it smoothly over the closed 
streamlines. Similar difficulties arise if a non-uniform shear flow is considered. 


x_x,b «a? 
2b°2C Dr 206% 


. 
1 


ae ae Q 2b? 


f) 0 a 
L DO2D ER 200 
Slse le 
Be BE: C=Vo2+b? 
& # > h 


Fig. 6. Vortex pair in shear flow near wall 


$ 5. No such simple treatment is possible when there is a pair of vortices at the same distance 
from the wall. Each separately would require a different inclination of the shear-flow direction. 
If the distance 2 b between the vortices is large compared to a, the flow near either one is pre- 
sumably of the type just indicated. Accordingly a pair which produce upwash between them- 
selves (fig. 6) would converge in the shear-flow direction. 

On the other hand, experimental evidence [2] shows that a and b are comparable and a > b 
in certain eircumstances, so that the mutual interaction of the vortices cannot be neglected. When 
»/a2 k is small, see (8), it is reasonable to assume that this interaction may be computed as if the 
vortices remained parallel to each other and the wall. Fig. 6(b) shows the velocities induced at 
the vortex 1 by the vortex 2 and the image vortices 3 and 4 under this assumption. The compo- 
nents are always positive as drawn. The only other velocity with which the vortex 1 tends to 

11 


- 
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' i i the vortex must therefore 
ove is that of the main shear flow, effectively®). In order to be free t L N 
Gera slightly inwards and upwards from the shear-flow direction, the direction cosines being 


au ” x a2 
Jaallı) Ib U(a)’ 


i i isi iri d [2], though no quan- 
spectively. This convergence and rising of the vortex pair is also observed [2], tho 

titative een are available. A further conjecture is that their disruptive action will continually 
increase as they rise. = 1 

The asian of groups of longitudinal vortices may be conjectured similarly. In particular, 

an infinite array of vortices, equally spaced at the same distance from the wall, wi of equal but 

alternately opposite rotation, will tend to converge in pairs as before but not rise). - 
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ara ae  Ee 
Kia ee a? +(2 m b)e 


= 1 x 


nı@®+(@n—ı0® #5 


2 a 
—ıa (est 7@ _tanh 3) e 


a 
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Der Spannungszustand in einer geschichteten Scheibe”) 
Von Hans BUFLER 


Es wird das zweidimensionale, ebene Problem der Spannungsbestimmung in einem aus einzelnen homo- 
genen und isotrop-elastischen Schichten mit beliebigen Elastizitätskonstanten und beliebigen Höhen bestehenden 
Körper bei gegebener Randbelastung oder Randverschiebung mit Hilfe der Fourier-Transformation all- 
gemein gelöst, und zwar sowohl für den Grenzfall vollkommener Haftung als auch den Grenzfall verschwindender 
Reibung zwischen den Schichten. Als Beispiele werden die homogen-geschichtete Halbscheibe und die Halb- 
scheibe mit anschließendem Streifen behandelt. 


By means of Fourier transformation a general solution is found for the plane problem of determining the 
stresses in a body consisting of homogeneous elastically isotropic layers of arbitrary elastic moduli and arbitrary 
heights when the boundary load or the boundary displacement is given; the solution applies both to the limiting 
case of perfect adhesion and to the limiting case of zero friction between the layers. Examples treated are 
the homogeneously laminated half-plane and the half-plane with adjacent strip. 


Ilpı moMmomm mpeo6öpasorannsı Dypne pemtaerca B OÖMeM Bunte ABYXMepHaf TNIOCKAaA 
3aylada OlPETeJIEHHA HANPSDREHMÄ B TEE, COCTOAMIEM NZ OTIIENIBHBIX OAHOPOAHEIX H H30TPOHHO 
DIACTHUHBIX CJIOEB € IPOUSBOJIBHLIMH MOAyıAMa JOHra u UPONU3BONBHON BbICOTbI PH 3anamoüi 
HArPy3Ke MM IIPU 3AaHHOM CMeINeHHUN Kpas, A HMECHHO KaK MIA CAIYyyaa HOAHOTO NPHJIHM- 
HaHnA, Tak M UI CJIyYası HMYTOKHOTO TPeHHA M&RAY OTACNBHLIMH CAIOAMM. B KayecetBe 


UPHMePOB PAccMaTpHBalmTcaH OJNHOPOAHO CIONcTag IOJIYITIACTHHRA M TNOJYNJAACTUHRA C 
IPUMBIRAIOINeH K HEH TIONOCON. - 


i 


1. Problemstellung 

Der Gegenstand der vorliegenden Untersuchung ist die 

in einer aus einzelnen Schichten aufgebauten e 
nichtperiodischer) Randbelastung bzw. R 


Bestimmung des Spannungszustandes 
lastischen Scheibe bei gegebener (periodischer oder 
andverschiebung. Als Schichten werden dabei lauter | 


*) Von der Fakultät für Maschinenwesen und Elektrotechnik der Technischen Hochschul ü 
genehmigte Habilitationsschrift (Referenten: Prof. Dr. L, Förpt, und Prof. Dr. H. No) Erin sa ee 
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unendlich lange, homogene und isotrope Streifen mit beliebigen Elastizitätskonstanten und be- 

liebigen Höhen, jedoch gleichen Breiten zugrunde gelegt. Ist die Breite der Schichten sehr viel 

größer als ihre jeweilige Höhe, so wird sich im wesentlichen ein ebener Formänderungszustand 
einstellen; ist die Breite dagegen klein im Vergleich zur Schichthöhe, dann kann man von einem 

ebenen Spannungszustand sprechen. 


Die zwischen den einzelnen Schichten auftretenden, vom jeweiligen Reibungskoeffizienten 
abhängigen Reibungskräfte haben maßgebenden Einfluß auf die Größe und die Verteilung der 
Spannungen. Es lassen sich dabei zwei Extremfälle unterscheiden: Der Grenzfall vollkommener 
Haftung und der Grenzfall verschwindender Reibung, bei dem nur senkrecht zu den Schichten 
wirkende Kräfte möglich sind. 


Nach Aufstellung der Lösung für einen beliebig belasteten einzelnen Streifen mit Hilfe der 
Methode der FouRIEr-Transformation werden die Berühr- und Haftbedingungen (letztere ent- 
fallen bei Reibungsfreiheit) formuliert. Diese ergeben sich als System linearer homogener Gleichun- 
gen für die Transformierten der Randspannungen der Schichten. Für einige Sonderfälle werden 
geschlossene Lösungen angegeben. ’ P 


Die Behandlung des rotationssymmetrischen Belastungsfalles eines geschichteten Körpers 
ist auf analoge Weise mit Hilfe der HAnker-Transformation möglich und soll einer anderen Arbeit 
vorbehalten bleiben. 


2. Mathematische Vorbemerkungen 


Bei der Untersuchung des Spannungszustandes in einem unendlich langen Streifen mit 
parallelen Rändern (weiterhin kurz Streifen genannt) kommen — je nachdem eine periodische 
oder nichtperiodische Randbelastung vorliegt — FourIzr-Reihen oder FOURIER-Integrale vor. 


Für die FOuURIErR-Reihe der periodischen Funktion /(x) mit der Periodenlänge 2 ! gilt — falls 
1 


+ 
/(«) im Intervall (— 1, I) die DiricaLer-Bedingungen!?) erfüllt und [ f(x) dx absolut konvergiert — 
die komplexe Darstellung?) 1 


mit 
0] Ku).e- ° "du ER a IE 
Nach Einsetzen von c._„in die erste Gleichung erhält man 


+! 


Q Be. u—x 
2 2. [we Kan dus) (1). 


an 


en 
tur) 1 
—i 


1 ELIA : 
An einer Sprungstelle liefert (1*) statt /(x) den Wert D [/(& + 0) + fx — 0)]. Die (eindimensio- 


nale) endliche FourıEr-Transformierte der Funktion f(u) (oder I&)) sei in ihrer komplexen Form 
als von der diskret Veränderlichen n abhängige Funktion definiert durch 


Kd=37, 
ei] 


+1 


mn TN 


ER - 
Fr EOS | fir) -e de 3 aa oe. 


7 
Gilt nun für f(x) die Darstellung (1*), dann folgt in Verbindung mit (2*) für die inverse Trans- 


formierte der Funktion f(x n/l) 
+ 00 


1 - {zn 1 2 P 


= — 0 


\ ; ee : Er X En SERRT. 
S Funktion f(x) erfüllt im Intervall (a, b) die DirıicuLrr-Bedingungen, wenn siein (a, b) eine beg 
eo Merums 2 Mlinims und eine begrenzte Anzahl von endlichen Sprungstellen — dagegen keine un- 
endliche Sprungstelle — besitzt. 
2) Vgl. z.B. G. P. Toustow, Fourierreihen, Berlin 1955, VEB Deutscher Verlag der Wissenschaften, S. 33. 
LVZ 


160  H. BurLer, Der Spannungszustand in einer geschichteten Scheibe 


und, je nachdem f(rn/l) eine gerade (d.h. f(an/l) = j—-zn/l)) oder eine ungerade Funktion 
(d.h. (anf) = — anf) ist, 


410 + S 17)" = Ir) gerade 


_1/2 a ac or 
W441 
re En ee - sin ” E (7) ungerade 


Schreibt man in (1%) A—= .n/l und führt anschließend den Grenzübergang 1— oo.durch, 
so geht diese Gleichung formal über in die komplexe Form der Fourıerschen Integralformel?) 


© + 
1 ENT 4* 
= —. DB 4 FF OU, 0er SE 
(d- 5, [a [io e d a*) 
welche sich für die Darstellung der nicht-periodischen Funktion f(x) eignet. Voraussetzung ist 
+ 
auch hier, daß f(x) für — © <x< + oodie DiricHter-Bedingungen erfüllt und f /(@) dx ab- 


— 

solut konvergiert. Ferner liefert (4*) an einer Sprungstelle statt f(x) den Wert 1/2 [f(x + 0) 

fe — 0)]. Definitionsgemäß ist nun die (eindimensionale) Fourter-Transformierte der Funk- 
tion f(u) (oder f(x)) in komplexer Schreibweise 


+8 oo 
iM) = er | No) TE - | AH). Maag RL (5*). 


Im Gegensatz zu n handelt es sich bei A um eine kontinuierliche Veränderliche. Falls für dieFunk- 
tion f(x) die FOURIERsche Integralformel gültig ist, erhält man daraus in Verbindung mit (5*) die 


inverse FOURIER-Transformierte der Funktion f(A): 


2 © 


. | (Op ertas da: RE ... (6*). 


u 


1 
I@) = az 


Je nachdem, ob f(4) eine gerade oder eine ungerade Funktion ist, ergibt sich aus (6*) 


De N 


Du Sta) cosAx- aA f(4) gerade 
fe) V: | 
de Se if) -sinAx-dA (A) ungerade 
0 


| Die FOURIER-Transformierte der Ableitung einer Funktion läßt sich — wie man durch 
partielles Integrieren der rechten Seite von (2*) bzw. (5*) zeigen kann?) — unter der Voraus- 
setzung, daß die ersten (r —1) Ableitungen von /(x) für x — + !bzw. für — + ooverschwinden, 


darstellen durch 
le &e En n au) 
[re m 


bzw. 
| AN = (Tara ji 
Dabei bedeuten 
H 
(mn | . d’ f(x) le 
EEE STE! . 
/ =) 7 | Ta NNe dB in 
und in 
1 oo 
Ne | lee > 
(r) u SEE NET r 
9) 75 ee 1. 


?) Wegen des Beweises s. z.B. 1. 


Mc. Graw-Hill Book Comp. ING, S. 15. 
4) Ebendort, 8. 27. 


N. Sneppon, Fouriertransforms, New York/Toronto/London 1951, 


| 
| 
| 
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3. Spannungen und Formänderungen eines Streifens 


Als Voraufgabe für die elastizitätstheoretische Behandlung eines geschichteten Mediums 
muß der Spannungs- und Verformungszustand eines einzelnen Streifens, dessen Randspannungen 
vorgegeben sind, bestimmt werden. 

Es gelten folgende Spannungs-Dehnungs-Beziehungen: 


ee A (a), 
=, -5 u (a a) Tr RE Ban (1b), 
ee Me ee 
= N BO ER?. (id), 
WEgtn gm. VE NE ET . (le), 
en = dene N RES ERE RL 


x, y, z sind die kartesischen Koordinaten, u, v, w die Verschiebungen in x-, y-, z-Richtung, o,, 
0, 0, die Normal- und 7,,, Tyz, T;, die Schubspannungen. E bedeutet den Elastizitätsmodul und v 
die Querdehnzahl. Liegt ein ebener Formänderungszustand vor, so reduzieren sich die Gl. (1a) bis 
(1f) wegen & = yys = Yır = 0 auf 


ou 1 


ame Po, OT re ER 
gell ee on 
et Rh ee 
Beim ebenen Spannungszustand erhält man dagegen wegen 09, = 7, = T,,; = 0 
nl Ba ee 
Spare ER N 


Der Vergleich der Formeln (2) und (3) zeigt, daß man die entsprechenden Beziehungen des 
ebenen Spannungszustandes aus denen des ebenen Formänderungszustandes dadurch gewinnt, 
daß man (bei festgehaltenem E) die quadratischen Glieder von » streicht. Hätte man statt 
des Elastizitätsmoduls E den Schubmodul G gemäß G = E/[2 (1 + »)] eingeführt, so käme man 
von den Formeln des ebenen Formänderungszustandes zu denjenigen des ebenen Spannungs- 
zustandes dadurch, daß man (bei festgehaltenem G) statt » den Ausdruck »/(l + ») setzte. Wir 
legen im folgenden die Gleichungen des ebenen Formänderungszustandes nach (2) zugrunde und 
können die dafür abzuleitenden Formeln mit »® = 0 sofort in diejenigen des ebenen Spannungs- 
zustandes überführen. 

Zu den Formänderungsbeziehungen (2) bzw. (3) kommen noch die Kräfte-Gleichgewichts- 
bedingungen am Element d - dy hinzu. Diese lauten (sowohl für den ebenen Formänderungs- 
zustand als den ebenen Spannungszustand) bei Abwesenheit von Massenkräften: 


00. OT 
02 0y \ (da), 
ee ee 


162 H. BurLer, Der Spannungszustand in einer geschichteten Scheibe 


Hierbei wurde 77, = Tys — T gesetzt, was aus der Momenten-Gleichgewichtsbedingung am Ele- 


ment folgt. r 

Der Streifen erstrecke sich über den Bereich — o<x2< + ©, 0<y<sh (Bild). Stellt 

man die Spannungen und Verschiebungen (bei nicht-periodischer Belastung) durch ein FOURIER- 
sches Integral entsprechend (6*) 


| MD = nz | ay-etrd ... 0) 


X 
Bild1. Der unendliche Streifen mit parallelen Rändern 
dar, so lassen sich die Randbedingungen längs y = 0 und längs y = h befriedigen. Die in (5) 
als Parameter erscheinende Größe y muß deswegen eingeführt werden, weil die Spannungen und 
Verschiebungen sowohl von der Koordinate x als auch von der Koordinate y abhängen. Die 
FOURIER-Transformierte zu f(x, y) lautet gemäß (5*) 

+8 


Ken a Te N 2 
y2r 


— 09 - 
Im Fall einer periodischen Belastung sind statt (5) und (6) die Transformationsformeln 
(8. (8*) und (2*)) 


n 


= u 5) 
und 
ne ji r“ 7 
f 7’ 779 I&,y)-e dx (6’) 
7 


zu verwenden. 

Die Spannungs- und Verschiebungskomponenten müssen (beim ebenen Formänderungs- 
zustand) die partiellen Differentialgleichungen (2a), (2b), (2c) und (4a), (4b) erfüllen. Multipli- 
ziert man diese Gleichungen mit e'** und integriert über die gesamte x-Achse, so erhält man 
nach (6) bei Beachtung von (8*) für die Fourrer-Transformierten folgende fünf gewöhnliche, 
lineare und homogene Differentialgleichungen mit konstanten Koeffizienten?): 


1 


MA y)= nr MN) - OMAN] - - - - (M 

2 1 ei | 
Dv(),yJ)= 5 (1 —9)0,(A,))— +) o,(l,Y)] - - - : &, 
Düß, ia) = any). re ee 
—Ae(A,)HDUÜD=O . 0 
DoylAW.— TAT). 


ee et den Operator d/dy. Genau auf dieselben Gleichungen kommt man bei Ver- 
x ; = ransformation (6'), sofern man in (7) bis (11) A durch z n/l ersetzt. 
Fre de Se Regen: können vorstehende Gleichungen nur dann be- 
ee ah var izienten bzw. Operatoren der Spannungs- und Verschiebungskom- 
En Ben aebagle erminante, angewandt auf irgendeine der Spannungs- oder Verschiebungs- 
groben, verschwindet. Diese Bedingung führt zu 


er (D — 220, J)=0. Er RN. 

V ] NA -\ ay% pe Te Er ee 7 

ee Kuh EnEwedez In Oyı T, U oder v gesetzt werden kann. Die Lösung der Differential- 
8 (12) lautet mit G,, C,, C, C, bzw. A, B, C, Dals „Integrationskonstanten‘“®) 


9% y) = G(A) - eflv+ C,): e-Wlv + C,(@): Al y- erlv + Ca) - |aly-eWlv, 


5 V |. ji ier; 3 T \ n < Y . 
Verlg.s. 74. ee BON SR D.S. Berry im Handbuch der Physik, Bd. 6, Berlin 1958, Springer- 
Minaer Methode ER unendliche Halbscheibe bei nicht-periodischer senkrechter Randbel st Fr h 
’ , bei der sich die Einführung einer Spannungsfunktion erübrigt, behandelt mr 
’ 2 2 7 


) / 3 a 1} L ß diese Größen unabhä neie von der RK oordina- 
W 1r schreiben „Konst ınten {-} 

- D 1 en Re ‚un damit anzudeuten, da D> 

te Y sın ’ dagege n sind Ssıe ım allgemeinen eıne Funktion von j | | 
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oder in anderer Form geschrieben”) 
9@, y) = A(A) - Coj [A| y+ BA) - Sin A| y + C@) [al y-CojjAly-+ DA): Aly-Sinlaly (13). 
Setzen wir jetzt o, = g, so erhalten wir der Reihe nach aus (11), (10), (7) und (9) 


Be Me U a A ET ALIEN (14), 
Ay) =— . AN. 2.0 a a Be N (15), 
0,(), y) = — — a Wa. Are er a a N RE (16), 
ie i D: e 

u) =— 3 ß ++ | ee (17), 
v(},y) = ne +7 — DE — (1— 2) "| Ten (18). 


Die in (13) auftretenden ‚‚Integrationskonstanten‘“ A, B, C, D müssen durch die Rand- 
bedingungen festgelegt werden. Bei gegebener Randbelastung handelt es sich um das sog. erste 
Randwertproblem; als bekannt vorausgesetzt sind hier längs y = 0 die Druckverteilung p(x) 
und die Schubverteilung g(z) und längs y = h die Druckverteilung p(x) und die Schubverteilung 
q(z) (Bild 2). Die Randbedingungen lauten also: 


Dix) 

12,0) = do); t@,h) = —g@); lllie 

0,2, 0) = — P@&); o,(&, h) = — p(@); YAPOABNFE 

oder transformiert E h ea 

Et. gen, (19a), 1% 
WO) = EN... (19b), on 1} Fr 
er ae (19e), 2 | 
6, ( E h) ER D ( N) Ba ( 19d)- Bild 2. Randbelastung des Streifens 


Verwendet man nun g nach (13) in (14) und (15), so führen vorstehende Randbedingungen zu: 


ee Eon, en N (20a), 
a a a ER TARERARER - Urs (20b), 
ae — ME AS+BC+CMIRS+Q+D(MnC+S] .. . (2) 
Due rare + ehIC+ Due]: Tn.. 2 en (204). 


Zur Abkürzung wurde dabei 
SinAlh=&; Cfjlhr=E 
geschrieben. Nach Elimination von A und Bin (20c) und (20d) mittels (20a) und (20b) folgt: 


— (Ah+SOp+(AlhC+S)p igei-i ni Alh- A 
= Sim ..(20e), 
& B-MRSP-iy, es SOG +inARe 94 
n= Al" ER! ee 
© — (Ah)? 


”) Da — mit Ausnahme der ‚‚Integrationskonstanten“ A, B, C, D —in (13) ) stets eine positive Größe 
bedeutet, während die sonst vorkommenden / sowohl positiv as auch negativ sein können, sind dort an den 
entsprechenden Stellen (auch da, wo es nicht unbedingt erforderlich wäre, z. B. bei Coj |A| y) Absolutstriche 
angebracht. 
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Wir setzen nun die gefundenen Ausdrücke für A, B, C, Din (13) und dann (13) in die 
Formeln (14) bis (18) ein. Sodann erhalten wir für-die transformierten Spannungs- und Verschie- 
bungskomponenten, wenn wir zur Abkürzung 

j 


FerT 


einführen: 


5,4, y) = an ((&2— (a HB) p&ofaly— Lian + SHOP —(ahE+S)Pp ——- 


+ j@Ahd4 +j|Alh Sg] Sin A| y 
+ [da] + SO P—(ARC+S)B+jS*4+ Ah Sg) A yEoflAly 
+ [EB + alnSp-+ (ah —SOF— iA] h C-S)qllaly Sinjalyy - - - U), 


& 1 s er 4 u 
#9) = Sg AN iNInSp+ MR SYi-MhE SA 


x Sin |aly + (&® — (A.1)2) q Cof Al y 

+ [-jdalh+SOB+i(AhC+S)P+&4+]2h Sg] al y Sin || y 

+5 ihn Sp + MR SO B— (AnC—S)allaucof alu} (OR 
= — |] = = 7 
N) = + a2 hlnep—2ilaln— 94 

+2; (2]h& — &) a] Cof Al y 

+ [- (A| h+ SC) B +(A|hC+E&)B+j (an —2&)7—j]2]hSg] Sin [A| y 

+ MR + SHOP +MRE+ SP —-i&4—jalRSg] alyEoilaly 

+ [EB — ASP — lan SOg+ilalnc— SallalySinlalyt. . - (3), 


a 1 e 4 

EAU) = gl KA—r— 29) & + (+) Bd +27 —")AhSp 
— 21-9) (Ah=-SOd +21 — (ARE — S)7] Eoj |A!y 
+ li —9— 2%) (Ah + SO P— ill —»— 2m) (A RC+S)P 
+A+N AM — 21 )SHYG— (1 —r— 2) ST] Sin [Al y 
++ SO P—ilAlrc + SP S&g—AlnSz]ı + ]AlyCofAly 
+1 &B + j | Sp — (Ah SH, g+(A|nC-- Sg] +) ]AlySinalyt (24), 

ee 1 k Ri 

E.ı-v(A, y) = Sam If (1—9) (Ah+ SO P—2U1—W)(ARC+S)P 


+ KU 2) + +NAmyg+2j(— Ah Sq]Cof | y 

+ KU + am — 21-9) SHE +l—v— 2) AhSp 

Hit 92) (AR -SOg— j(1—v— 292) (A| AE— S)g] Sin [Al y 

+ AnSp— ja SOE+ RC Sal ıtn)]alyCofaly 

| + (Ah +SOp + (A RC+S)B— 1 &27— j|A| ST] I-+v) A| y Sin alt (25). 
Die beiden letzten Formeln gelten für den ebenen Formänderungszustand. Mit »? = (0 gehen 


sie nach dem eingan =) G Ag ) 1 


OA 

Di Baar die 1 eh (6) zugrunde, die bei nicht-periodischer Be- 
: :ine periodische Belastun , ir ie T i 

Sue a ee stung vor, so wird man die Transfo j 

en, welche ebenfalls zu den Formeln (21) bis (25) führt, sofern man dort die ee 


Veränderliche A durch die diskrete Veränderliche ersetzt 
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Die transformierten Randbelastungen p, 9, p, g, mü | } 
r 5 9, P, q, müssen aus den gegebenen Funk 
p, q gemäß (5*) bzw. (2*) errechnet werden: E BEBENEREN UNS lIoneN Dr 


DA) +0 [P@&) 

N ou 

Ba) = Yan Bi ae a N a (26) 

10) ie) 

(bei nicht-periodischer Randbelastung) 

) (p® 

| 1 leo |. 

en Sr | Bun ET ee (26’) 
an —ı P(&) 

(r) ia) 


(bei periodischer Randbelastung) 


Durch Rücktransformation gemäß (5) bzw. (5') erhält man 


0,(%, Y) [ (4, y) 
OR, Y) ı . 7 [@&»9 
ta, y), = Vor } | ATS) nee (27) 
ur, n —o | u(2, 4) 
v(X, y) v(A, y) 
(bei nicht-periodischer Randbelastung) 
ne | 
0(2 9) er 
{ren 
oy(z, y) 0, T » v) 
1 is REN =, 
rg yaı : 3 er - „) a A (27) 
NnN=—w 
u, y) r, „) 
v(x, y) Zur ; v) 


(bei periodischer Randbelastung; hier ist in (21) bis (25)4 = T zu setzen.) 


Damit ist der Spannungs- und Verformungszustand des pi) 
Streifens bei gegebener normaler und tangentialer Rand- || | | | | | | H 
belastung allgemein bestimmt. Lt uer Balz 


Bevor wir uns der geschichteten Scheibe zuwenden, 
wollen wir noch die Spannungskomponenten der homo- 
genen Halbscheibe aus der Lösung für den Streifen ab- 
leiten. Führt man in (21), (22), (23) die Koordinate n ge- 


mäß Y = h NN ein (vgl. Bild 2) und geht dann zur Grenze Bild 3. Die Halbscheibe bei normaler und 
h — oo über (Bild 3), so erhält man für die Halbscheibe tangentialer Randbelastung 
An) = — [+ AMP) — in aaler"?, 
an) =— ia Amar", 2.2... (8) 


(47) = — [ — lm) PA) — 7a — Am aAler "" 


Pan 


4% | $ ee 
2) = — RE [Ir — Alm) joa N) lea 


Für i im Intervall — a<x< a gleichmäßig verteilte HandepeniipE po und % folgt nach ( 


F En m 2 sinAa ee Ei S: 
ge: gr a A et EB Do: Mn eo 

E Po) y2r | N; K ge E L ce 

E | er er 
Ber ur 2 Sl Ad 3 "285 oe . 
Be. ? w=\2 Io’ 7 i E e . . 2.70 


x Geht man schließlich zu einer Schneidenlast mit den Komponenten P, = ar Ally 2a] und 


— 
(= lg [95 : 2 a] über?), dann erhält man nach a des Grenzibergasgs a—0— 
” a Be 0 „ 
] P= —- en m nr 30). 
p Vom q Vom (30) 


Wir setzen nun diese Ausdrücke in (28°) ein und bekommen bei Beachtung von (6’*) und Ver- 
wendung der Integrale 


je} 3 
2 
[ermrans fon PTR 


—1a| A ER dena 
ER "csArd= 2 ae rs DE 2 


en (31) 
e-lln gi wi Ag; BRRT.. 
sinAzdA=2|eA"snAzxd =, 
2 d 
+ x oo 
— A £ 23 —ÄAY Ü 9 xın 
r. InsinAzd=2|%e a re 
— 0 0) 
die bekannten Formeln 
u RR n 7? 2 2m? 
EDS ET re 
2 x? zen 
zu — I ® mm nn nn 
tx, n) - P, En = Dee @ En ei ER (32). 
, Ba a an 2 a8 
ine Po ar nip — 


+ 
9) P, bzw. Qyi 


Bi ist die im Koordinatenursprung angreifende Normal- bzw. Tangentialkraft je Scheiben- 


j 
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4. Die gesehichtete Scheibe im Falle vollkommener Haftung 


Wir betrachten jetzt einen aus N einzelnen Streifen bestehenden Körper gemäß Bild 4 und 
setzen voraus, daß keine Relativverschiebung der Schichten längs der Berührungsflächen möglich 
ist. Dieses trifft sicher dann zu, wenn die Schichten zusammengeklebt sind. Falls sie frei aufein- 
‚ anderliegen, soll der jeweilige Haftreibungskoeffizient so groß sein, daß keine tangentiale Ver- 
schiebung eintreten kann; ferner soll in diesem Fall eine geeignete äußere Lastverteilung das 
gegenseitige Abheben der Schichten verhindern. Die Schichthöhen A, h,,... Ay seien beliebig 
groß, ebenfalls die entsprechenden Elastizitätsmoduln E,, E,...Ey und Querdehnzahlen ER 

Ygy...vy. Vorgegeben seien die Randlasten, 
nämlich p,(x) und q,(&) an der oberen und py(«) 
und qy() an der unteren Begrenzung (s. Bild 4). 
Gesucht sind zunächst die ‘zwischen den 
einzelnen Schichten übertragenen Normal- 
und Schubspannungen :Pı(®), 4ı(&) Ps(x), 
IX), ::-» Pr-ı&), gv—ı(2). Die Vorzeichen 
von pund gsind dabei durch Bild 2 festgelegt. 


Drrs, Pk Pr-ı, 
Ir «I 


Bild 4. Die aus N Streifen zusammengesetzte Scheibe Bild 5. Zur Herleitung der Berühr- und Haftbedingung 


Zur Herleitung der Berühr- und Haftbedingung betrachten wir den k-ten und (k + 1)-ten 
Streifen (Bild 5). Die Bedingung dafür, daß sich diese beiden Streifen nicht voneinander abheben, 
lautet offenbar: 


v,(&, 0) = Yp+ı (&, Ar+ı)  (Berührbedingung). 


Soll ferner keine tangentiale Verschiebung der Berührungsflächen auftreten, so muß außerdem die 
Bedingung E 


1,(2, 0) = Ug+ı (%, A,4ı) (Haftbedingung) 
erfüllt sein. Wir gehen jetzt zu den transformierten Größen gemäß (6) über und erhalten?) 


DKM ON= BET (A, dur) ° (Berührbedingung) . 7 3. 22°.7...@83), 
HA, O) 223. (A, Ripr).# (Haftbedinsung) Tessa a Ar He 


Für die in (33) und (34) benötigten Verschiebungen ergibt sich nach (24) und (25), wenn man dort 
PD; DB» q; g durch die entsprechenden Ausdrücke Pı-1, Pr Pr+1 Ik—1, Ik Qr--ı ersetzt: 


ler). 


11.(4, 0) = zer MDR t+tik ra — Entf at - -:-*. 85), 
IpeneH 5 2 
Ari „2 Ir = A a“ £ a 4 
>. a0 > = ai 1% RR — Br ty et Ti}. 0. 302286), 
Exp & 
2(1 — +1) : = m FB *IRRT 37) 
= —— + {7 Öp+ı Pr+ı — ] Ya+ı Pe — Bir Ierı + Ofrı A} (57), 


Alu) 
A» Eryı * &+1 


Der (Ay hr) = - 1Pr+ı Par — +1 Pr + TÖrtı erı + I YıHılı 8). 


°) Die Transformation (6) ist bei nicht-periodischer Belastung, die Transformation (6') bei periodischer 
Belastung anzuwenden. Im zweiten Fall braucht A lediglich durch ı n/l ersetzt zu werden, <o daß sich eine 
gesonderte Behandlung dieses Falles erübrigt. 
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Dabei wurden folgende Abkürzungen verwendet: 
&% = %(4) = A], + Sin A| hr  Eof |A| Ir, 

art = arlA) = |A| hr — Sin A| hr - Cof j2| Ar, 

Br = Pr) = |A| Ar - Cof |A| x + Sin [A| Ar, 

Bk = Pia) = |A| Ar - Co] [A| hy — Sin jA| hr, 

dr = Ör(A) = |A| hr - Sin |A| Ir, 9). 

= — Sin} — (A hy)? 
f no %—2 SC + (i +) (A hy)? = (1 — 2») Sin?Ahy + (A hr)? 


1— 1—% 
ae beim eb. Formänd. Zustand 
(1—») Ein A ha + (1 -+ 94) (A iu)? beim eb. Spannungszust. 


Mit (35) bis (38) gehen die Berühr- und die Haftbedingung über in 
Br+ı * Pexzı — (r+ı + 0% * Pr, n+1) Pr + Pr Pu,x+ı * Pr-ı 


+ jögrı -Aerı + 1 Ya — Pr Pt) m —T Pr ka=0 . . . (0), 
(Berührbedingung) 
Ör+ı * Pr+ı — (Yr+ı — Ye Pr,k+1) "Pk — Op "(Pk,k+1" Pr-ı 
+ iBkrı > Terı — I orı + 0ER - PunrHı) et IßE Pr ma =0 .. . (Al). 
(Haftbedingung) 


In vorstehenden Formeln bedeutet 


1—% Ertı &+ı 


beim eb. Formänd. Zust. 
a 


kt = Ip u PR TREE 
N paal beim eb. Spannungszust. 

BO & 

Die Gln. (40) und (41) sind zwei gekoppelte, lineare und homogene Differenzengleichungen 

zweiter Ordnung für die Funktionen p;(A) und 4,(4). Bei N Schichten sind 2 (N — 1) dieser Funk- 

tionen, nämlich p,, Pa - - - Pv—15 I» 9» - - - 1x—ı zu bestimmen (P, 9, und px, qy Seien durch die 

Randbedingungen festgelegt). Dazu stehen uns die GlIn. (40) und (41) zuf Verfügung, welche 

(N — 1)mal angeschrieben werden können (nämlich fürk = 1,2,... N — 1) und somit genau die 

erforderlichen 2(N — 1) Gleichungen liefern. Damit ist das Problem im Prinzip gelöst; um die 

Spannungs- und Verschiebungskomponenten an einer beliebigen Stelle der k-ten Schicht zu er- 
mitteln, hat man wie folgt vorzugehen: 


1. Aus den gegebenen Randbelastungen py(&), Qo(%), Pr(X), Ax(&) sind die transformierten - 
Belastungen p,(A), go(A), Pr(A), Qv(A) zu berechnen, und zwar bei periodischer Belastung nach 
(26°) (A = sn/l) und bei nicht-periodischer Belastung nach (26). 

2. Die 2(N — 1) Gleichungen (40) und (41) (für k=1,2,...N — 1) sind unter Berück- 


sichtigung der Randbedingungen nach den P,, Pa - - - Pv—ı5 Qi» des : » : Iv_—ı aufzulösen (im periodi- 
schen Fall ist dort A = a n/l zu setzen). 


5 
1 
- 


3. Nach Einsetzen von Pr, 9% Pr, Ik in (21) bis (25) erhält man die transformierten 
Spannungs- und Verschiebungsgrößen 6,, 0, T, u, v der k-ten Schicht (wobei man im periodischen 
Fall = zn/l zu schreiben hat). 


4. Durch Rücktransformation gemäß (27) (im periodischen Fall) bzw. (27) (im nicht- 
periodischen Fall) gewinnt man die Formeln für die Spannungen o,, 0, t und für die Verschie- 
bungen u, v an einer beliebigen Stelle der k-ten Schicht in Summen- bzw. Integraldarstellung. 


Falls am oberen und unteren Rand der geschichteten Scheibe statt der Spannungen die 
Verschiebungen (zweites Randwertproblem) oder teils Spannungen, teils Verschiebungen (ge- 
mischtes Randwertproblem; die Spannungen oder Verschiebungen sollen jedoch hier jeweils im 
gesamten Bereich — 00< x < oo gegeben sein) vorgegeben sind, hat man in dem für die erste 
Schicht (k = 1) und die letzte Schicht (k = N) geltenden Bedingungsgleichungen (40), (41) die 
Transformierten der Randspannungen p,, G; Px» 4x ganz oder teilweise durch die Transformierten 
der Randverschiebungen 1,(A, h,) = u,(4), D; (A, h,) = vu(A) und uy (4,0) = Uuy(A), Dy (A, 0) = Dy(A) 
auszudrücken. Dies geschieht für die erste Schicht mittels der für k — 0 angeschriebenen Gin. (37), 
(38), für die letzte Schicht mittels der für k — N angeschriebenen Gln. (35), (36). Ist beispielsweise 
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am oberen Rand die Verschiebung v,(x) vorgegeben, während dort keine Schubspannung auf- 
treten soll (9, = 0), so folgt aus (38) mit k = 0 und D, (A, h,) = v,(A) sowie 9, = 0 
= 1 E,e& = & 

> 4 At. 


Damit lauten die Bedingungen (40), (41) fürk = 1 


en 2 ei Eee 4 0) = 
Pa P2 -(# +0 nen) 2/0 ment ea) 
1 


Pı &ı E, ee f 
cu 91,2 Da a NER Ar a at, Mon ee (40’); 
= Dr N ER BR " ö7 2 
6, PL — [Ye — Yı Pu. a Pı +1 a i(e + 6 "Paz Pa di 
1 
8 E R 


2 ee, 


a 2m) 
Wir betrachten im folgenden einige Spezialfälle. 


4.1. Die homogene Halbscheibe 

Die Halbscheibe bestehe aus lauter gleichen, aneinandergeklebten Streifen. Die Lösung 
der Differenzengleichungen (40), (41) muß dann auf die bekannten Formeln der Halbscheibe 
führen. Dieses wollen wir im folgenden nachweisen. 


Wegen der Gleichheit der Schichten wird nach (42) or,2+1 = 1; außerdem ergibt sich nach (39) 
ee hieß eier mi VE Ma=ot: Be Drsr= B* 
Die beiden Differenzengleichungen (40), (41) erhalten mithin konstante Koeffizienten: 

PPıs1 — 20 Pr + PP 10 > On De, (43) 
a dra tin lernt. nn. ii). 
Zur Lösung verwenden wir die Ansätze 
BesATet ug esBNer 
welche folgende zwei lineare, homogene Gleichungen für die „Konstanten“ A und B ergeben 1°): 
A-(B-e&e —2ae+Pß)+B-öleer —1)=0 
A.s(&e— 1) + B(f* &e —2ate + Br) =0. 
Eine nichttriviale Lösung existiert nur dann, wenn die von den Koeffizienten gebildete Deter- 
minante verschwindet. Diese Bedingung liefert die charakteristische Gleichung für o und lautet 
mul ee 2: 
(BB*— 8) 4-2 (aß* + Bar)? +2(BB* + 200 +0) 2 — 
— 2( $ß* + Bar) z + —-P—=0. 
Setzt man jetzt die entsprechenden Ausdrücke nach (39) ein, so folgt nach kurzer Umformung 
(2 —2Coj Alh-z2 +1? = 0 
mit den Wurzeln , = 3=ett;, =2,= el, (entsprechend =, =lhh a =4,=—J[Ah). 
Damit haben wir die Lösungen 
Dr(A) = Ay(A) - eMlhkE Aa). Alhk. ellkL Ad). ei 4 A,()-|Ahk.erite, 
144) 32,0). re 4 BO): Ahk-efl®® BA). e-@l»rı BO).|Aalhk-e-iihe, 
Um den Zusammenhang zwischen den B und den A zu finden, setzen wir (45) in (43) ein und 
nehmen einen Koeffizientenvergleich vor: 
(BCof A Ak— a)-A, + hBSinlA|hk- Aa+6&in|Alhk-B, + A hökof[Alhk-B,=0, 
(B&of A hk—a)-A,+66&inAlhk-B,=0, 
(BCof AlAk— a): A, — A| hßSinAjhk- A, —6Sin|Alhk. BB+|A hößojAlhk-B,=0, 
(BCojAl Ak— a) -A,—65&inlhk.-B,=0. 


| is 


10) Die Bezeichnung ‚‚Konstanten“ soll andeuten, daß diese Größen unabhängig von k sind, während sie 
im allgemeinen von A abhängen. 
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Mit Berücksichtigung von (39) erhält man schließlich 2 
Br De N 
B, FE A, 5 B, = A, - 

Die „Konstanten’’ A,, Ay, Ag, A, werden durch die Randbedingungen festgelegt. Da für k — oo 

die Spannungen verschwinden müssen, folgt aus (45) A, = A, —=0. Für k = 0 kommt dann aus 

(45) mit (46) 

Po) =A, und (4) = A; — A, ; 
Bei vorgegebenen Randspannungen (erstes Rand- 


wertproblem) sind p,(4) und qg,(4) als bekannt zu 
betrachten; mithin ergibt sich für die „‚Konstanten‘“ 


A, und A, 
A;=Pp,(A) und A, = Ppo(A) — 7 g0(A), 


Bild 6. Die Halbscheibe mit aufgeklebtem Streifen . 
belönorngpleXBelsstung und wir erhalten aus (45) 


DA) = [1 + Ah) po) — I Al hkgo@)]  eTRH 2 
4) = I |Alhkpua) + (1—|A| Ak) go(a)] - e="* nn (AM. 


- Diese Ausdrücke sind identisch mit den in (28) für die homogene Halbscheibe direkt abgeleiteten 
Formeln, sofern man dort o, = — pP, T= — q, undn = hk setzt, was nachzuweisen war. 


4.2. Die Halbscheibe mit aufgeklebtem Streifen 


Wir betrachten den aus einem Streifen (Index 1) und aus einer Halbscheibe (Index 2) 
bestehenden Verbundkörper gemäß Bild 6 und setzen der Einfachheit halber einen ebenen 
Spannungszustand voraus. Es sei eine beliebige Normalbelastung p,(x) des Randes vorgegeben 
(von einer tangentialen Randbelastung wollen wir hier absehen). 

Mit A, — ooerhält man aus (39) für die Halbscheibe 


e RE RR CE e2 Al he KRERLER elAlhz | 
Ba nen BE R=-P=-%=|Alh, ER, 
e2 All 
»=(l— 9) Lage RR (beim ebenen Spannungszustand). 
Nach (42) wird dann 
E, el 


(beim ebenen Spannungszustand). 


Pr = —— u 
P1,2 Ede 


N > © 


Setzt man diese Ausdrücke in die für k = 1 angeschriebenen Differenzengleichungen (40) und 
(41) ein und geht dann zur Grenze h, — ooüber, so ergibt sich wegen 9, = 0 (da keine Tangential- 
belastung des Randes vorhanden sein soll) 


Da 1— Ey uB 
I 2a). a in S,.2, 2 ern 
| Er Al 2 * Be I E & Po » 


lu, Ey a ER Sn. 
| ee 


1 


Die Auflösung dieser beiden Gleichungen nach pı und q, führt zu 


E \ IR, “ 
© Pı° = be Be *) + du | 5 a = ae „) 
Pı(A) = PA): a a . 


E E ln Br 
1 =] nm \E ” 
re en a KA Bee 

= De +). € “) | 5) E, & n) 


1l—», E E 
: 2 Re ( ee on)tafgats 
m 3 ea 


2 5, ee 
wata)lga ) | >) na n) 


s 


Ber 
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Po(A) ist. aufGrund der vorgegebenen Randbelastung als bekannt anzusehen; die &,, af, e, und Yyı 
können aus (39) entnommen werden. Durch Rücktransformation von (48) und (49) erhält man 
die zwischen dem Streifen und der Halbscheibe übertragenen Normal- und Schubspannungen. 


Bestehen Streifen und Halbscheibe aus demselben Material, so vereinfachen sich (48) und 
(49) (der Nenner wird (Sin |A| h, + Eof |A| hy)? - [Sin2A h, — (A h,)?]) zu 
PıQ) = Po(A) (1 + Ahr) e=Mit, 
HA) = —iPo(a) - [A| hy, e= Mm. 
Diese Formeln stimmen — wie es auch sein muß — mit jenen der homogenen Halbscheibe über- 
ein. Davon kann man sich durch einen Vergleich mit (28) oder (47) überzeugen. 


Wir wollen noch einen zweiten Spezialfall ins Auge fassen. Setzt man in (48) und (49) 
E,— © (was dem Fall eines elastischen Streifens auf starrer, rauher Unterlage entspricht), so 
folgt 

+6 
Dee ‚Aa +4 Yı 
Pı( ) Po(A) %ı 0° + y? ’ 


— Pıyı FH 
My +Yi 


In Verbindung mit (39) erhält man daraus nach Division von Zähler und Nenner mit 
[Sin®ih,— (A hy)?] 


(4) = TPo@) - 


& En 2CojAhh+(l+v)AhSinAh 
pi — AN: 1 3 1 1 
PıA) Po(A) 40 Al, — (1 —v) Sin®Ahh+(1l+ 9m)? (Ah) 
r „= (1 —») Sinkh—(l+»v)Ahß&ojiAh, 
A) 2 . 3 
aa ZEHN an Ad Emil ÜTMÜhR 
Falls nur eine (normal zum Rand gerichtete) Einzelkraft P, bei x = 0 wirkt, gilt nach (30) 
= 1 
— hi Ze P , 
Po y2 m 0 
und die Rücktransformation von p, und gq, gemäß (6’*) liefert 


ep 0. 260fAh+(+v)Ah&inih, 
AMT 70° | dCi A PEmihtlitr®An) 


0 


cosAx:dA (50), 


oo 


_ 51 — $ 
eg = Pe | (lv) Snih—(+v)Ahbojih, ee, (51). 


460)? 7 h, — (1 — Sin’}hh+ (1 +9)’ (Ah) 
v0 
Dieser Fall des elastischen Streifens auf starrer, rauher Unterlage unter Einzellast wurde 
früher schon von K. MARGUERRE !!) behandelt. Seine Ergebnisse stimmen mit den Formeln (50) 
und (51) vollständig überein. 


5. Die geschichtete Scheibe im Falle verschwindender Reibung 


Wir legen unseren Überlegungen wieder einen aus N Schichten von beliebigen elastischen 
Eigenschaften bestehenden Körper zugrunde, wie ihn Bild 4 zeigt. Der Unterschied gegenüber 
dem unter der Ziff. 4 behandelten Fall liegt darin, daß keine Reibung zwischen den einzelnen 
Streifen auftreten soll. Mithin ist nur eine senkrecht auf.die Schichten wirkende Belastung mög- 
lich. Diese soll so verteilt sein, daß sich die Streifen nicht gegenseitig voneinander abheben. 


Die Bedingung für Reibungsfreiheit lautet 


Geh elr=I N ET A 

Sie ersetzt die Haftbedingung, die hier entfällt. Die Berührbedingung muß dagegen erfüllt 
werden. In Verbindung mit (52) reduziert sich die Berührbedingung (40) auf 

Br +1 Przı —(a+ı + % Omar ı) Pr + fr Oma+ı  r-ı=0... . (8). 

Gl. (53) ist eine lineare und homogene Differenzengleichung zweiter Ordnung für die Funk- 

tion P,(A). Da insgesamt N Streifen vorhanden sind, müssen (N — 1) unbekannte Funktionen 


11) K. MARGUERRE, Ing. Arch. 2 (1931), S. 108. 
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Dis Pas - - - Pv_—ı bestimmt werden (p, und py seien durch die Randbedingungen festgelegt). Zu 
diesem Zweck steht uns die Differenzengleichung (53) zur Verfügung, welche, fürk=1,2,... 
N — 1 angeschrieben, gerade die erforderliche Zahl von N — 1 Gleichungen liefert, womit die 
Aufgabe im Prinzip gelöst ist. Hat man auf diese Weise z. B. p; und p,—ı bestimmt, so erhält 
man aus (21) bis (25) die transformierten Spannungs- und Verschiebungsgrößen 9,, 0,, T, u, v der 
k-ten Schicht und nach Rücktransformation schließlich o,, o,, T, u, v. 


Ist statt der Druckverteilung p,(x) die Verschiebung v,(x) vorgegeben, so gilt für k = 1 die 
Berührbedingung (40) mt = 9% =: 
a. | - 2 € ‚E m ie 
Bed (mt ma fa) mar ag An .. . BI). 
Wir wenden uns jetzt einigen Sonderfällen zu. 


5.1. Die homogen-geschichtete Scheibe 
Sämtliche Streifen sollen hier gleiche elastische Eigenschaften und gleiche Abmessungen 
(Höhen) besitzen. Dann gilt 
= +4ı =; er g=&%41=8&; 9,r+ı =1 
und (53) geht über in folgende Differenzengleichung mit konstanten Koeffizienten: 


Per 2 Pr + ni =0 2 Re 


Der Lösungsansatz 
Pr EN 
führt zu 


& 
ee —2-e+1=0. 
ß 


Daraus folgt 


& [a\2 | 

edı,2 — + IN 1 4 

PT VB ö 

bzw. ’ 

& (&\2 ’ 

a=—a=nls + /6)-ı)- j 

| P ß ö 

Nach (39) ist aber s 

& _Ah+ &inAhlofah _ j 

B Ahcofiht Smih= 

und somit ist 9, 20. Die allgemeine Lösung lautet nun mit A,(A) und A,(4) als an die Rand- 

bedingungen anzupassende „Konstanten“ ü 

Bud) = Ad) - at + Ayla). mat, 

oder, nach Einsetzen von p, und r 
u) 


wobei PN) AHA nennen. 
yo) = Ah+ SinAh-Cojah+ySin®A h — (Ah)? - SinAh 
An-Cofäht&mih 0 0 
bedeutet. 
Mitk=0undk=N lautet (55): 


Po = A, | A, S Px : A, R yN + A, .p7 N ; 


x pP nd l V SEeDe ‚ Bi 5 Ss oO S R l a zusehen. Die Auflösune 
ar 8 n 


A, == Px a Me “Po . A ER YN Do — Pr 
uN — werde Sr ie yN— wN . 
und damit wird 
Pr(A) En ne al N N > k . en = 
| gr Bw... 0. 
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Der Übergang zur homogen-geschichteten unendlichen Halbscheibe ist schnell 


vollzogen. Setzt man nämlich in (57) N — oo, so folgt sofort 
DEAD A a 


Wir geben jetzt noch die sich durch Einsetzen von y aus (56) und durch Rücktransformation ge- 
mäß (6*) ergebende Formel für die Druckverteilung zwischen der k-ten und der (k + 1)-ten 


Schicht an 2): 


00 


Ah-CojAh+ SinAh “eg (58). 


e. 1 2 Ah+&inAh-CofA h+yS©in®A h— (A h)2-SinAh]l * 
Pr) = Yan Po) - 


—09 


Für eine Schneidenlast P, am Rand ist in bekannter Weise p, — P,/Y2 r zu setzen; allerdings 
müssen wir in diesem Fall noch eine — wenn auch geringe — gleichmäßige Druckvorspannung 


überlagern, welche ein evtl. Abheben der Schichten voneinander verhindert. 


FürAh<1 geht (56), wie man durch Reihenentwicklung nach Potenzen von Ah zeigen 


kann, über in 


(An): 
| B-1+,m 
und (58) wird in diesem Fall 
(Ah)? 
£ = Ah 21—%k 5 ne 
Pr(A) & Po(A)- 1" En en 3 Die Alten a ue 


der Form 


P+ı—-2r+m-ı 2/0 TEN 
h2 n2\ 8 1):r=0, 
wenn man zur Grenze A—0, d.h. zur Differentialgleichung 


dp [2/a Er 
et 2 


übergeht. Durch Reihenentwicklung nach Potenzen von A h wird 2 en — 


Diesen Ausdruck bekommt man übrigens auch aus der Differenzengleichung (54), geschrieben in 


und die Lösung der vorstehenden Differentialgleichung liefert bei Berücksichtigung der Rand- 


bedingungen 
nn _ am: n 


a 
und mit n = h-k sofort den Ausdruck (59). 


5.1.1. Erstes Beispiel 


Als Beispiel bestimmen wir im folgenden den Spannungszustand in der homogen-ge- 


schichteten Halbscheibe für die Belastung (vgl. Bild 7) 


Bla) p*.: ( + cos a 1.60) 


Nach (26°) erhalten wir für die transformierte 
Belastung 


Tr 
V5 . p* Tre + 1 
Ds = für (607). 
y2 1a p* n=0 


Bild 7. Die homogen-geschichtete Halbscheibe bei periodischer, 


0 sonst normaler Belastung 


Damit wird (58) = 
l 
V» k 
0 


: | 
r)- ST: pr Be Kin nn 


0 ( sonst 


12) Bei periodischer Belastung ist natürlich mit! = n/l die Transformation (3*) zu verwenden. 
12 


TE 


nr 
4 


r 
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n 
wobei nach (56) wegen eh 


ar, Sin + ein Ver #-) 


3 


7 
„[e® en 
l sch sach UNSER 
| Te „62. 
n 1 /ahı2 Eh b 
er Be). 
yo) =1. 
Die Rücktransformation von (61) gemäß (3’*) liefert 
Die) — pr: ( + Ir 7) a ee ER (61). 


Für relativ dünne Schichten (d.h. die Schichtdicke h soll im Vergleich zur Perioden- 
länge 2 ! sehr klein sein) erhält man mithin 


ah\2 & 
1. fahhalm m | Se "ays 7 
SD 1 + — | — .cos—ıı vv p* il 17.605 — 61’). 
DAS PD hl] con ART ER cos 7x ( ) 
Der rechtsstehende Ausdruck in (61’') ist identisch mit dem von G. SonNTAG ®) angegebenen, 
der durch eine sehr vereinfachte, speziell auf den Fall dünner Schichten bei periodischer Rand- 
belastung zugeschnittenen Ableitung erhalten wurde. 

Wir wenden uns nun der Bestimmung des Spannungszustandes in der k-ten Schicht zu und 
wollen der Einfachheit halber z h/l < 1 voraussetzen. Nimmt man in (21), (22), (23) 


P=P; P=M-ı=v'Pr; i=d20, 
so folgen mit hl <1 


[+00] (2H)] er Se 

‚ RER 

n=6j eo. W y|! n): im): 2. (08 

3 eh 

ee o( 2) rn Dee NE 
”) 


und nach Rücktransformation entsprechend (3’*) und Einsetzen von Y 


Er 1 N I. aha 2 T, y 7 
nn] Banane retro en 


=. | 


2.2.64), 


f 
ı 
L; 
h 
i 


2 En 


022(%, y)= y3 ib Pi 


Bild 8. Verlauf der Spannungen über die Streifenhöhe in der homo 


gen-geschichteten Halbscheibe mit x h/t<1 (Belastung nach Bild 7) 


Bild 8 zeigt den Verlauf von Oyks 


7%; und 0,7; über die Streifenhöhe. Man erkenn rei i 
Biegespannung ist; der Verl ee 


auf von r;, entspricht ebenso der von der elementaren Balkentheorie 


13) G. SonNTag, Forsch. N Se EHE 
ots a ‚ Forsch. Ing. Wes. 23 (1957), 8. 3 ; dort werden die Schichten nach der elementaren Balken- 


J 
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her bekannten parabolischen Schubspannungsverteilung. Falls die hier vorausgesetzte Bedingung 
schjl< 1 nicht erfüllt ist, wird natürlich eine mehr oder weniger große Abweichung der Span- 


Bild 9. Verlauf der Spannungen in der homogen-geschichteten Halbscheibe und in der nicht-geschichteten Halbscheibe, 
aufgetragen über % (Belastung nach Bild 7; m h/l = 1/10) 


nungen von den in Bild 8 gezeigten Kurven vorhanden sein. Von Interesse sind noch die maxi- 
malen Spannungen im k-ten Streifen, aufgetragen über k. Aus (63°), (64°), (65) folgt: 


2I—k+1 
oa fit ]143,(®)] hrs) 
arsch |, 1 Fahren 
use tale] .p*, EEE TECHN 


ee 1 sch 2T—k 
te] 


Zum Vergleich geben wir noch die leicht aus (28) ableitbaren Formeln für die Spannungen in der 
nicht-geschichteten Halbscheibe an. 


We 1 + v 4% 2.) Eu .p*, 


(Dep =T ne on”, ee nr iv): 


Ber 
(O7)2=0 — | -- ( 7 ) 6 s | . DR 


Bild 9 zeigt den nach (66) und (67) für z h/l = 1/10 errechneten Verlauf der Spannungen, auf- 
getragen über k (in (67) wurde dabei n = kh gesetzt). Die für (66) zuständigen Kurven geben 
natürlich nur an den diskreten Stellen k=1,2, 3, ... die richtigen Spannungen (und zwar die 
maximalen Werte für die k-te Schicht) an. Man erkennt das im Vergleich zur nicht-geschichteten 
Scheibe außerordentlich langsame Abklingen der Randspannungen in der geschichteten Scheibe. 


5.1.2. Zweites Beispiel 

Hier soll der Spannungszustand in der durch eine Schneidenlast P, am Rand beanspruchten 
homogen-geschichteten Halbscheibe untersucht werden. Für den Druck am unteren Rand der 
k-ten Schicht gilt (s. (58°)) 


je) 


Biz)— = Iv® BEINE A, ee N) 
Fe 
mit 


1* + Sin A* » Co]2* + Sin? A* —1*?-©ini* 
uch Ar CofA* + SinA* ; 


176 H. BurLErR, Der Spannungszustand in einer geschichteten Beheibe 


wobei 4* und x* die dimensionslosen Größen A* =Ah und x* = z/h bedeuten. Die tangential 
gerichtete Normalspannung am unteren Rand der k-ten Schicht beträgt 

u © 5 

P | 2a y- Smart —(SindAr FA), u. cost zr-d* . . (69). 


Rn Denn Sinzir — m 


y . 
Diesen Ausdruck erhält man durch Rücktransformation von (23), wenn man dort y = 0 und 


g= {= 0 sowie D = Pı_ı und P = Pu setzt: 


1 
Sin?A* — Ar? 


5,1(4*, 0) = {24% : Sin A* - Pa_ı — (Sin?A* +4*9). Di} - - (70). 


h 


0.102: 3 4 3, I EDEN 


— 


Bild 11. Verhältnis der Hauptschubspannung in der geschich- 

teten Halbscheibe zur Hauptschubspannung in der nicht- 

geschichteten Halbscheibe längs der Wirkungslinie der normal 
gerichteten Schneidenlast (gültig für k=1, 2, 3,...) 


Bild 10. Verlauf der Spannungen 27.(0), 0,7(0,0) und 
Ty (0,0) in der durch eine normal gerichtete Schneidenlast be- 
anspruchten homogen-geschichteten Halbscheibe längs der Wir- 
kungslinie der Kraft (gültig für k=1,2,3,...) 
er Verlauf der entsprechenden Spannungen in der nicht- 
geschichteten Halbscheibe 


Bild'12.,— Verlauf der Spannungen 9,(2*) und 0% 1@*, 0) in 

der durch eine normal gerichtete Schneidenlast beanspruchten homo- 

gen-geschichteten Halbscheibe längs des unteren Randes der 1. Schicht 

- Verlauf der entsprechenden Spannungen in der nicht-ge- 
schichteten Halbscheibe 


an rn ri LE EEE GNS 


Die Integrale (68) und (69) lassen sich in geschlossener Form nicht streng lösen. Sie wurden : 
für x* = 0 (d.h. längs der Wirkungslinie der Kraft P,) numerisch ausgewertet; das Ergebnis . 
j ! ze r 1 
zeigen die ausgezogenen Kurven von Bild 10, wo auch die Hauptschubspannung 7; = p2 (Pr + 0zx) i 
eingetragen ist. Zum Vergleich sind die entsprechenden für die nicht-geschichtete Halbscheibe | 
maßgebenden Kurven (gestrichelt) eingezeichnet, für welche | 
er ee 
Pie (0): =0 = Py: n ’ (G.-o=0, (n = kh) 


gilt. Man erkennt — wie beim ersten Beispiel — das im Vergleich zur nicht-geschichteten Scheibe 
ziemlich langsame Abklingen des Drucks infolge der Schichtung. Bild 11 zeigt das Verhältnis der | 
Hauptschubspannung in der geschichteten Scheibe zur Hauptschubspannung in der nicht-ge- 
schichteten Scheibe (gültig jeweils für den unteren Rand der k-ten Schicht und längs der Wir- 


% 


H. Burter, Der Spannungszustand in einer geschichteten Scheibe 127 


kungslinie der Kraft P,). Sieht man die Hauptschubspannung als Maß für die Anstrengung des 
Materials an, so folgt, daß die Schichten im Vergleich zu den entsprechenden Stellen in der nicht- 

F geschichteten Scheibe um so stärker beansprucht werden, je tiefer sie liegen. 
Ferner wurde für den unteren Rand der am stärksten beanspruchten ersten Schicht (k=71) 
der Verlauf von p,(x*) und o,1(2*, 0) aus (68) und (69) ermittelt und in Bild 12 niedergelegt, wo 


a die Spannungen für die nicht-geschichtete Halbscheibe zu finden sind. Letztere errechnen 
sich aus 


2 Bo a Ne A En 
mh (ae oa SOmsaTe nn Gare 


Man erkennt insbesondere die sehr hohe tangentiale Randspannung o,,. 


p(«*) => — (den = 


Bild 13. Spannungen p,(O)undo,, (0, 0) in dem durch eine Schneidenlast normal belasteten, auf einer Halbscheibe 
reibungsfrei liegenden Streifen in Abhängigkeit vom Parameter f 


5.2. Die Halbscheibe mit anschließendem Streifen 


Dieser Fall unterscheidet sich von dem in Ziff. 4.2. behandelten dadurch, daß keine Schub- 
spannung zwischen Halbscheibe und Streifen wirkt. 
Mit h, — 00 bekommt man aus (39) für die Halbscheibe 


e2lAlhz elAlhz 
a nina, I Ps = A| hz 2 |, 
a > © n a > © 


Damit wird aus (42) 


1-2 E, ellh 


Zr (beim eb. Formänderungszustand), 
912 = 1 = 


—» E 4a In-o 
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und die Differenzengleichung (53), angeschrieben für k = 1 und aufgelöst nach Pı, geht über in 


Pl) = Bo 7 E 
DE ER 


Mit (39) folgt daraus 
% # A| hy - Cof || A, + Sin [A| h, j 
De | 3 
2 Er. (Sin? A, — (Ah)? + Ah + Sin JA] Ar «of A] 
2 


2 
1—v] 


Besteht die Belastung aus einer Schneidenlast P, bei =, so gilt nach (30) D = P,V2r a 
und die Rücktransformation von p, gemäß (6’*) liefert unter Voraussetzung eines ebenen Form- 
änderungszustandes (beim ebenen Spannungszustand braucht man nur »—=»3=(0 zu setzen) 
folgende Druckverteilung zwischen Streifen und Halbscheibe: 


oo 


P, 7* .CofA* + Sin A* De ee, ER 
Pıl@*) = on, | 7.7 Simeir —AM) FAR + Sin Ar. Copar °0° 71) 
0 
1-2 E 
Dabei bedeuten A* = Ah, a* = x/h und [= —— 2.1. 


(Eine gleichmäßig verteilte Zusatzbelastung des Streifens soll auch hier dafür sorgen, daß sich 
dieser nicht von der Halbscheibe abhebt). Setzt man in vorstehender Formel f—= 0, was dem 
Fall eines elastischen Streifens auf starrer, glatter Unterlage entspricht, dann bekommt man 
Übereinstimmung mit dem von E. MELAN'4) schon früher erhaltenen Ergebnis. 


Bild 14. 


Verlauf der Spannungen p,(x*) und Oy,(a*, 0) längs des unteren Randes des durch eine Schneidenlast normal belasteten, 
auf einer Halbscheibe reibungsfrei liegenden Streifens für /=1 
ee Verlauf der entsprechenden Spannungen in der Halbscheibe ohne Schnitt 


Für die tangentiale Randspannung am unteren Rand des Streifens folgt nach (70) mit k = 1 


oo 


pP 1 
ar, Va = 
age sch, ler Ar — 1*2 
0 
ee (Sin? A* + 4*2) (A* CofA* + SinA*) 
x 124% Ol A# „Ar FREENET ee . cosA* x*. dA* (72). 
| F- (Sin2A® — Ar) + A* + SinA®.Gopanj SoSAttdA* (72) 
Die numerische Auswertung der Gl. (71) und (72) zeigen die Bilder 13 und 14. In Bild 13 
ist die Abhängigkeit des unter dem Streifen auftretenden Maximaldrucks Pı(0) und der maximalen 
Spannung 0,,(0, 0) vom Parameter f dargestellt. Bild 14 gibt den für /=1 (d.h. Streifen und 
Halbscheibe bestehen aus gleichem Material) errechneten Verlauf von Pıl<*). und o,, (2*, 0); 


zum Vergleich sind wiederum die entsprechenden Spannungen in der homogenen Halbscheibe 
eingetragen ®). 


14) E. MeLAn, Beton und Eisen 18 (1919), 8. 83. 
'5) Siehe hierzu auch K. MARGVERRE, Ing. Arch. 4 (1933), 8. 332. 


j 


a 


Schichten. 
weise 
Ajı Aıa 
1 Aı Agg Ags 
Ag Agz Agy 
Ags Ay 


wobei 


Ar,k—ı = Pr* Ok,k+1 5 
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6. Anhang 


Zur Lösung des Gleichungssystems (40), (41) und des Gleichungssystems (53) bei beliebigen 


Das für den reibungslosen Fall zuständige Gleichungssystem (53) lautet in Matrixschreib- 


Pı a, 

P3 0 

Bar er 
Ay; 0 ; 


\ An-ı 


Ap,c+1 — Pr+ı 


An—ı,n—2 Ay-ı, ne oe 


Ak = — (r+ı + Or * Pr,0+1) ; 


und bei vorgegebenen Randlasten 


a, 


Pr’ ua Po; an ı = — Pu "PN 


(in den freigelassenen Feldern stehen Nullen). 
Es läßt sich durch fortgesetztes Eliminieren (GAussscher Algorithmus) auf folgende Form bringen 


b,ı bis Pı b, 
baa ba; Pa b, 
Da Da, Ps en b; 
bja Dis Pa b, 
ee) Px-ı bu) 
mit 
a : er Ik—1,5 ' Ag,r—ı . 
Dr — djı P ber = Ark z b 3 bk—1,% — dk—1,k (k == 2 3 ah N ee 1) 
k—1,k—1 
=; Belt... as (k=2,3,...N—2) 
b,ı Da2 bk-1, k—1 
dd AN—1.N—: 
br = NW at. 22... 2 + ayoı. 
b,ı Da5 Dn_2, N—2 
Daraus erhält man der Reihe nach 
r by = r 
Bra =, — 5 Bra g—— (nr — dran noanı Pre) (K=23,...N—D). 
N-1,N—1 N—k, N—k 


Für das Gleichungssystem (40), (41), das für den Fall vollkommenen Haftens gilt, wird man 
keine geschlossene Lösung anstreben, da hierbei außerordentlich umfangreiche Ausdrücke ent- 
stehen. Vielmehr wird man das Gleichungssystem (40), (41), das von der Form 


(Ay, A152 d13 1a 
Ay Aga Ayz Agı 
Ay] Ayy Ayz Aza Ag; Age 
Ayı Ayo Ay Aya Ays5 Ass 
Agz Azy Az; Age Asy7 Ass 
Agz Agg Ag5 Age Ag7 Ags 
Ay; Age A77 Ay A79 Ar, 10 


ER: 
W, Ag 
w,; 0 
w, 0 
W, ie : 
Ws; äh 
w, 
0 
daN—3,2N—5 QaN—3,2N—4 dan—3,2N—3 d2aN—3,2N—2 | | Wan—-3 d2N—3 
dan—2,2N—5 d2aN—-2,2N—4 A2N-2,2N—3 ee) W2N—2 Aa EG 
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ist, für eine Reihe von 4- 
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Werten nach den je 2 (N — 1) Unbekannten ;(A) [pı(4) = W2r-1(A); 


j gu(A) = War) (k= 1,2, ... (N — 1))] mit Hilfe des Gaussschen Algorithmus numerisch auf- 


lösen 1%). 


Abschließend sei bemerkt, daß das Problem der geschichteten Scheibe sich auch mittels 


Übertragungsmatrizen behandeln läßt, worauf an anderer Stelle noch eingegangen werden soll. 


16) R. ZurmÜühL, Prakt. Mathematik, 2. Aufl. Berlin 1957, Springer-Verlag, S. 91ff. 
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KLEINE MITTEILUNGEN 


Heat-Treatment of a Sphere made of Visco- 
Elastic-Plastic Material with Temperature-De- 
pendent Yield Stress 


We consider a sphere made of visco-elastic-plastic 
material being heated or cooled down from a known 
temperature T = Tr). The temperature distribution 
at any subsequent moment is assumed to be known. 
The problem is to find or, o, and 09 as functions of 


r and £. We shall neglect the effect of stresses on 
temperature distribution and the inertia forces. 

For the visco-elastie flow only, we write the stress- 
strain relations as 


; 5 1 Ichivı. : 
=) +) WM. 
Er 1 2 € Or + 2 er) 
3 =. 3 TEEN N (2): 
The equation of equilibrium is 
Be 
ne =. 0 fein Me (3) 


and the condition of compatibility 
0 
a) Dr .(4). 


We have assumed in the above that the tempera- 
ture distribution is a function of r only sothat e, = &9, 
H=pandoyg=0,,= == ee. 
For the portion of the sphere wherein the plastie flow 
is not zero, equations (2), (3) and (4) hold good, while 
instead of (1) we have 

(#r — 05)? = Y&T) 
so that in this region 


YH—-,=+]7 


: 5) 


. (6) 


and 
4 —u=+Y ne Rear 


At any given moment, a given point of the sphere 
is eitherin a plastie stage (&; = 0) orin a visco-elatsie 
stage. If it isin a visco-elastic stage, it stays in this 
stage so long as |o| < Y. We want to decide under 
what conditions does a change of phase, from visco- 
elastic to visco-elastic-plastic or vice-versa, take 
place. 

So long as a point is in’ visco-elastie stage, equations 
(1), (2), (3) and (4) apply and we can caleulate sand & 
from these equations. If, at the moment these equa- 
tionsgive o = Y, theyalso give 6 > Y, the point will 
undergo a change of phase from visco-elastic to visco- 
elastic-plastie and from that moment onwards, till 


the point stays in this latter stage, equations (2), (3), 
(4), (6) and (7) will apply. 

Similarly, if at the moment equations (1), (2), (3) 
and (4) give a=—Y, they also give 6<—Y, 
the point will undergo a change of phase from visco- 
elatic to visco-elastic-plastic and from that moment 
onwards, till the point stays in this latter stage, 
equations (2), (3), (4), (6) and (7) will apply. 

The conditions under which a change of phase will 
take place from visco-elastie-plastic to visco-elastie 
are as follows: 

Suppose the plastic flow started at a given point 
at a moment t=t,. Also let us suppose that at this 
moment i=t,, at the point under consideration o is 
equal to Y, as given by equations (1), (2), (3) and (4). 
Then these equations also give at this point, at the 
moment t=t, 


> . (8). 
Plastie flow continues at this point for moments 


t > t, only so long as the value of ö given by (l), (2), 
(3) and (4) satisfies inequality (8). This means that 
the plastic flow takes place, only so long as ei; = (0 
will give 6> Y (at the points where «= Y). As 
soon &,= 0 will giveö< Y, the plastic flow will 
cease, 

Similarly, if plastic flow starts at a point at the 
moment t = t, and at this moment, at the point under 
consideration, equations (1), (2), (3) and (4) give 
6 = — Y then they also give at this point at t=t, 


en Re 


. As before, plastie flow will continue at the point 
under consideration for moments t > t,, only so long 
as the value of ö as given by equations (1), (2), (3) and 
(4) for moments t > t, satisfies inequality (9). 


We shall see presently that these conditions satisfy 
the requirement that during plastic flow, the plastic 
work Se; >0. 

We notice that equations (2), (3) and (4) are valid 
at all moments at all points of the sphere. 

Substituting from equations (3) and (4) in equation 


(2) and integrating, and observing that &, is finite at 
r=(, we get 


Ä 
Kart |Brtedr 2 LO: 
0 


Substituting from (10) in (1) we get 


r 
1% 1 Da 
Hetze -ergale [onen) (11). 


0 
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Equation (11) gives for those points where |o| < Y 


Ü BL t 
= exp (- [e#) jr. t) = [? u)a (12), 
{0 


where 


er, =2r —. E | 3re “) MERLa). 


ö 
At those points where o= + Y, equation (11) 


gives 
m 
) 
We may now verify the condition of positive plastic 
work in the plastie regime. From (10) 


5 Pfr) + . (14). 
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giving, at those points where Ei; +0 
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0 
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accordingaso=+Y/. 
The condition of positive plastic work is 


&P __ ;P 
=> 0 


accordingaso=-+Y/. 
This condition reduces to 


1 
34 fi : : Y 3(1—») ; 
BEN 2 = 
„|?rear+serg, = H 7 
) 


which is equivalent to the condition that 6, as given by 


(14) should be greater than Y or less than — Y 
respectively. Hence, etc. 


Suppose we consider a sphere being cooled down 
from any temperature and we want the stresses in 
the sphere as a function of r and t. Let us suppose 
that tt=0,1 =,=0,=0. 

Let us fix our attention at any point r. In the be- 
ginning equation (12) will give the value of o at this 
point. Suppose that at #= t, equation (12) gives, at 
this point o = Y. 
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If at i = t,, equation (14) gives, at this same point 
6= Y, equation (12) will still apply for moments 
> - 

If, however, at t = t,, (14) givesö > Yatthe point 
under consideration, then for 3 >t>t, we shall 
have,o = Y,6= Y where t; is the largest value of t, 
such that fort3 > t> t, (l4)givsöo> Y. 


Att=t,o= Y and for moments t > tz, 6 must bo 
caleulated from (12) and we get 


t t t | 
ol Fa) fonnem[[Fa)arr 
2 


(15). 

Equation«(15) will give the value of o at the point 

under consideration, so long as this value satisfies 
|o| = Y and so on. 


If we know o| = — ” r (or) |, a simple integration 
will give o, and then equation (3) will give o,. 


Nomenclature 


r, 06,9 = spherical co-ordinates with origin at the 
center of the sphere 
t = time 
A = Temperature 
Er, Ep Ey &,9 Erg Ep E strain components 
Or, Op Ogp O9 Org dp E stress components 
E,v,n= known functions of r and 
= (1—2»)/E - K is assumed to be constant 


Y = the yield value of |or —o,|}; Y is a known 
function of T 

(3 = 00 

&; = plastic components of deviatoric strain rate 


Si; = deviatorie components of stress 

E = [ad 

& = co-efficient oflinear thermal expansion 

dot == denotes differentiation with respect to time 
ß = E/[6(1—v)] 


Reference 


H. G. LANDAT, J. H. WEINER and E. E. ZwICKY JR., Thermal Stresses 
in a Visco-elastic-plastic Plate with Temperature Dependent Yield 
Stress. J. Appl. Mech., Ser. E, Vol. 27, No. 2 (1960). 


Verfasser: B.D. AGGARWALA, Rensselaer Polytechnic 
Institute, Dept. of Mechanics, 
Troy (New York), USA 


BUCHBESPRECHUNGEN 


F. Winckel, Technik der Magnetspeicher. 
XVI + 6148. m. 325 Abb. Berlin/Göttingen/Heidel- 
berg 1960. Springer-Verlag. Preis geb. DM 72,— DM. 


Als Gemeinschaftsarbeit führender Fachleute auf 
den verschiedenen Teilgebieten der magnetischen 
Speichertechnik bietet das vorliegende Werk in seiner 
Gliederung Theorie— Konstruktion — Messung — Nor- 
mung sowohl dem neu an dieses Fachgebiet heran- 
tretenden Wissenschaftler und Techniker eine gute 
Einführung als auch dem Fachmann auf dem jeweili- 
gen Fachgebiet wertvolle Anregung. 

Hier seien nur die Beiträge herausgegriffen, die für 
die Technik der Rechenanlagen von Bedeutung sind. 

H. Biıruıng (München) betrachtet die Magnet- 
speicher bei elektronischen Rechenmaschi- 
nen. Bei der prinzipiellen Darisgung des Magnet- 
trommelsperichers werden die verschiedenen Schreib- 


verfahren, Kopfformen und Auswahleinrichtungen 
behandelt, der Magnetplattenspeicher wird erwähnt. 
Umlaufspeicher werden in Form von dynamischen 
Trommelspeichern beschrieben. Ebenso stehen Nickel- 
verzögerungslinien, Magnetband- und Magnetkern- 
matrixspeicher an. 

Der Digitale Speicher als Bausteinin Daten- 
verarbeitungsanlagen stellt besondere Anfor- 
derungen an die Organisation der Speichereinrichtung. 
R. PıtLoTy (Stuttgart) zeigt die verschiedenen Identi- 
fizierungsverfahren hinsichtlich der Ausnutzung des 
Speicherraumes und der Aufrufzeit. 

Meßverfahren an Speichern für kontinuierliche und 
diskrete Informationen werden von W. KALLENBACH 
(Braunschweig) unter Meßtechnik des Magnet- 
speichers besprochen. Auf Gleichlaufmessung, 
Spaltjustage und Frequenzgangmessung wird einge- 
gangen. 
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F. Kronss (Leverkusen) behandelt in seinem Bei- 
trag Theorie des Magnetspeichers knapp aber 
verständlich statische Speicher. Drei große Ab- 
schnitte befassen sich mit dynamischen Speichern und 
behandeln den Wiedergabevorgang, die magnetischen 
Eigenschaften der Magnetschicht und den Aufsprech- 
vorgang. / h 

Die Zahl der selbständigen Kapitel und ihrer Auto- 
ren läßt die wissenschaftlich-technische Breite und 
Bedeutung dieses Speicherverfahrens erkennen. Es 
ist dem Herausgeber in hervorragendem Maße ge- 
glückt, durch die Heranziehung der verschiedenen 
Autoren und der damit verbundenen unterschied- 
lichen Blickpunkte eine fruchtbare Vereinigung der 
aus Theorie und Praxis erwachsenden Problemstellun- 
gen zu erreichen. Es ist dem Buch, das einem be- 
stehenden Bedürfnis nachkommt, zu wünschen, daß 
es auch als Lehrbuch an Hochschulen Eingang findet, 


Jena W. KÄMmMERER 


H.L. Dryden, Th. von Kärmän, and 6. Kuerti, Ad- 
vances in Applied Mechanics, Vol. VI. X + 
294 S.m. 63 Abb. New York and London 1960. Aca- 
demic Press. Preis geb, $ 9,00. 


Der Band wird von zwei großen Arbeiten beherrscht. 
In der ersten bringen G. Lupwıc und M. Heır (Boun- 
dary-Layer Theory with Dissociation and Ionization) 
eine umfassende Darstellung über die bei Höchstge- 
schwindigkeiten notwendige Weiterentwicklung der 
klassischen Grenzschicht-Theorie. Die zwei Teile der 
Arbeit sind fast unabhängig voneinander lesbar. Der 
erste beschäftigt sich mit der kinetischen Theorie 
dissoziierter 2-atomiger Moleküle (Gase A,); unter 
Berücksichtigung der Transportgleichung für mole- 
kulare Figenschaften von Teilchen A, ergibt sich der 
Wirkungsquerschnitt beim Zusammenstoß eines Ga- 
ses A, mit Atomen; dieser Teil dürfte wohl nur dem 
theoretischen Physiker verständlich sein. Der 2. Teil 
vermittelt jedoch dem Aerodynamiker wichtige Er- 
kenntnisse; hier werden die laminaren Grenzschicht- 
gleichungen bei einem dissoziierenden Gase A, aufge- 
stellt und deren Lösungen im Fall der ebenen Platte 
und der Staupunktströmung angegeben. 

Die zweite umfassende Arbeit stammt von Kr. Os- 
WATITSCH (Similarity and Equivalence in Compressible 
Flow), die auch für den Berechnungsingenieur von 
größtem Interesse ist. Verf. geht davon aus, daß die 
linearisierten gasdynamischen Gleichungen sowohl 
bei Unter- wie auch bei Überschall in einfachem 
Zusammenhang mit denen der inkompressiblen 
Strömung stehen; hierbei ist bei der Transformation 
der Geschwindigkeitskomponenten ein konstanter 
Faktor A frei wählbar; systematisch werden Ähn- 
lichkeitsregeln aufgestellt, wobei affine Beziehungen 
zwischen Körper und Modell bestehen. Die Ergebnisse 
werden darauf ausgedehnt, daß auch quadratische 
Störungsglieder berücksichtigt werden. Bei den Äqui- 
valenzsätzen werden Körper geringer Flügelstreckung 
mit konischen Rotationskörpern verglichen. Die Ar- 
beit, in der Näherungslösungen mit exakten Lösungen 
und Meßergebnissen verglichen werden, schließt mit 
der bekannten Flächenregel. 

Die wertvollen Beiträge von K. SrrwArtson (The 
Theory of Unsteady Laminar Boundary Layers) und 
W. Chester (The Propagation of Shock Waves along 
Ducts of Varying Cross Section) zeigen, daß auf diesen 
Gebieten der Forschung noch ein weites Feld offen 
steht. Der Band schließt mit grundsätzlichen Bemer- 
kungen von R. Wırur (Kärmän Vortex Streets) über 
Wirbelstraßen-Probleme, wobei vor allem die bekann- 
ten Arbeiten aus dem Föttinger-Institut herangezogen 
werden. 


Eichwalde b. Berlin W. BADER 
K.P. Stanyukovich, Unsteady motion of Con- 


tinuous Media. XV + 7458. London/Oxford/Pa- 
ris/New York 1960. Pergamon Press. Preis £ 5 net. 


Das Werk behandelt eine Fülle von Problemen aus 
dem Gebiet der instationären Strömungen in gasför- 
migen Medien und der Ausbreitungsvorgänge von Stö- 
rungen in Wasser und festen Körpern. Mathematisch 
gesehen handelt es sich um Anfangs- und Randwert- 
probleme bei einem System von 2 oder 3 partiellen 
Differentialgleichungen 1. Ordnung; trotzdem sind die 
verwendeten mathematischen Hilfsmittel verhältnis- 
mäßig einfach, aber die straffe, souverän geschriebene 
Darbietung des umfangreichen Stoffes verlangt von 
dem Leser ein beträchtliches Maß von Mitarbeit. Das 
Werk stützt sich naturgemäß auf die diesbezügl. so- 
wjetischen Forschungen (hier sind außer K. P. Stanyu- 
KOVICH insbesondere L. D. LanpAv, L. I. Sepov und 
I. B. SELDoOVIcH zu nennen); es würdigt jedoch auch 
die bekannten Arbeiten von R.Covrant und K. 
FRIEDRICHS, K. BECHERT, E. JouvgvErT u.a. In dem 
Werk wird lediglich die analytische Lösung der Pro- 
bleme behandelt, wobei zwar die EuLersche Darstel- 
lung überwiegt, jedoch auch von der LaGrangeschen 
Darstellung in großem Ausmaß Gebrauch gemacht 
wird. Vorzugsweise wird von der allgemeinen Lösung 
— z.T. unter der Annahme recht allgemeiner Zu- 
standsgleichungen — ausgegangen. Das vom Aero- 
dynamiker bevorzugte Charakteristiken-Verfahren 
wird nicht verwendet. Die umfangreiche Monographie 
befaßt sich zunächst (200 Seiten) mit Problemen der 
eindimensionalen isentropen Gasströmung, wobei die 
Form k = (2r + 3)/(2r + 1) für den Adiabaten-Ex- 
ponenten angenommen ist. 300 Seiten etwa umfaßt 
die Darstellung der Theorie der stationären und in- 
stationären Verdichtungsstöße; auch die Ausbreitung 
von Detonationen wird berücksichtigt. Es folgt die 
Behandlung der dreidimensionalen Gasströmung und 
der instationären Bewegung in Wasser und festen Kör- 
pern; hierbei wird auch ein allgemeineres Fließgesetz 
den Rechnungen zu Grunde gelegt, so daß in diesem 
Fall die sonst unmöglichen Verdünnungsstöße auf- 
treten können. Anschließend werden (250 Seiten) Ein- 
zelfragen (Probleme der inneren Ballistik, Gasströ- 
mungen im Schwerefeld und bei verdünnten und sehr 
dichten Gasen sowie einige Fragen der relativistischen 
Mechanik) behandelt, wobei auch Fragen der Kosmo- 
gonie gestreift werden. — Dem Verlag gebührt größter 
Dank für die Herausgabe der Übersetzung. 


Eichwalde b. Berlin W. BADER 


Information Processing. Proceedings of the 
International Conference on Information Processing, 
Paris 15—20 June 1959. 5208. München 1960. 
Oldenbourg Verlag GmbH. Preis geb. 84,— DM. 


Es ist außerordentlich zu begrüßen, daß die zahl- 

reichen Vorträge, die auf der ersten Internationalen 
Tagung über Informationsverarbeitung (ICIP) ge- 
halten worden sind, mit den zugehörigen Diskussions- 
beiträgen und die in den einzelnen Symposien ge- 
lieferten Beiträge so schnell und in so vollständiger 
Form im Druck erschienen sind. Die UNESCO, auf 
deren Veranlassung und in deren Gebäude dieser 
weltweite Kongreß unter der Leitung von Prof. 
H. H. Arkev abgehalten worden ist, und die Verlage 
R. Oldenbourg (München) und Butterworths Scientific 
Publications (London) haben sich damit ein großes 
Verdienst erworben. Die Vorträge sind in mehrere 
Kapitel aufgeteilt, die den einzelnen Sektionen des 
Kongresses entsprechen und in ihrer englischen Be- 
zeichnung hier aufgeführt seien: 
1. Methods of Digital Computing; II. Common 
Symbolie Language for Computers; III. Automatic 
Translation of Languages; IV. Pattern Recognition 
and Machine Learning; V. Logical Design of Com- 
puters; VI. Special Session on Computer Techniques 
of the Future; VII. Miscallaneous Topics. 

Schon diese Aufzählung zeigt, — was durch die 
Themen der einzelnen Vorträge noch bestätigt wird —, 
daß sich der Kongreß eine sehr umfassende Aufgabe 
gestellt hatte und die damit zusammenhängenden 
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Probleme in ihrer ganzen Breite und Vielfalt erörtert 
worden sind. Auf Einzelheiten einzugehen verbietet 
die große Zahl der Vorträge. Es sei aber darauf hin- 
gewiesen, daß der vorgelegte Band durch ein ausführ- 
liches Stichwortverzeichnis (in englischer und fran- 
zösischer Sprache) und durch ein Verzeichnis der Ver- 
fasser zitierter Arbeiten zu einem willkommenen 
Nachschlagwerk wird, das in der Lage ist, jeden 
Interessenten über den Stand der Forschung auf 
seinem Gebiet zuverlässig zu informieren. Auch das 
Teilnehmerverzeichnis des Kongresses ist eine be- 
grüßenswerte Bereicherung des in jeder Hinsicht vor- 
züglich ausgestatteten Bandes. E 


Dresden H. HEInRIcH 


E. Schmidt, Einführung in die Technische 
Thermodynamik und in die Grundlagen der 


chemischen Thermodynamik. 8. Auflage. XX+ 
5438. m. 244 Abb., 69 Tabellen, 3 Tafeln. Berlin/ 


Göttingen/Heidelberg 1960. Springer-Verlag. Preis 


‚geb. 33,— DM. 


Das in Fachkreisen weit über die Grenzen Deutsch- 
lands hinaus bekannte ScHumiprtsche Lehrbuch zur 
Thermodynamik liegt bereits wieder in einer neuen, 
der nunmehr 8. Auflage vor, nachdem erst vor etwa 
2 Jahren, im Frühjahr 1958, die 7. Auflage erschienen 
war. Diese rasche Auflagenfolge (die 6. Auflage kam 
1956 heraus) ist nicht zuletzt ein Zeichen dafür, wel- 
cher großen Beliebtheit und Wertschätzung sich das 
SCHMIDTsche Lehrbuch in der Fachwelt, an den Hoch- 
schulen wie in der Praxis, nach wie vor erfreut. Der 
Unterzeichnete kann es sich daher ersparen, in dieser 
Besprechung auf die Qualitäten und Vorzüge dieses 
bekannten Standardwerkes, wie z. B. die konsequente 
Verwendung von Größengleichungen, die zahlreich 
beigefügten Tabellen, die in den Text eingestreuten 
und geschickt ausgewählten Übungsbeispiele usw., 
erneut ausführlich hinzuweisen. Das ist in den Be- 
sprechungen zu den früheren Auflagen bereits zur 
Genüge geschehen. Auch bezüglich der Aufmachung 
erübrigt sich jedes Wort. Sie liegt in bekannter Quali- 
tät vor. Vom Inhaltlichen her unterscheidet sich die 
8. Auflage von der vorhergehenden im wesentlichen 
durch die Hinzunahme eines Kapitels über den Zu- 
stand der Materie bei höchsten Temperaturen, den 
sogenannten Plasmazustand, sowie durch die Erwei- 
terung der im Anhang beigefügten Wasserdampf- 
tabellen auf den maximalen Parameterbereich von 
500 at und 800° C. Sie sind ein Auszug aus der vom 
gleichen Verfasser bearbeiteten und ebenfalls im 
Jahre 1960 erschienenen 5. Auflage der VDI Wasser- 
dampftafeln. 


Dresden N. ELSNER 

H. Schubert / R. Haussner / J. Erlebach, Vier- 
stellige Tafeln und Gegentafeln. (Sammlung 
Göschen, Band 81.) 1578. .Berlin 1960. Walter de 
Gruyter & Co. Preis brosch. 3,60 DM. 


Gegenüber der letzten Auflage von 1945 sind einige 


- Veränderungen vorgenommen worden: Fortgelassen 


wurden. alle ‚‚Gegentafeln“‘ (Antilogarithmen usw.), 
neu aufgenommen wurden Tafeln für Quadrat- und 
Kubikzahlen. So enthält das Bändchen jetzt folgende 
vierstelligen Tafeln: Dekadische Logarithmen der 
Zahlen 1 bis 10000 sowie der trigonometrischen 
Funktionen (Argumentschritt 1’ bzw. 10’), Additions- 
und Subtraktionslogarithmen, natürliche Logarithmen 
der Zahlen 1 bis 100, Quadrate und Kuben der Zahlen 
1,00 bis 9,99, natürliche Werte der trigonometrischen 
Funktionen (Argumentschritt 10’). Es folgen einige 
kleinere mathematische Tafeln (u. a. Umrechnung 
von Winkelmaßen, Zinsrechnung) und Tafeln mit 
geographischen, astronomischen und physikalischen 
Angaben. Hervorzuheben ist die übersichtliche An- 
ordnung, die durch schwarzgedruckte Numeri und 
rotgedruckte Logarithmen noch unterstützt wird. 
Dresden A. SCHUBERT 


0. Perron, Irrationalzahlen. 4. Aufl. VIIL + 
202 8. Berlin 1960. Walter de Gruyter & Co. Preis 
geb. 28,— DM. \ 


Das bekannte Buch des Verfassers hat seit seinem 


ersten Erscheinen im Jahr 1921 nichts an Frische und. 


Lebendigkeit eingebüßt. Daher ist auch die vierte 


Auflage gegenüber den früheren, besonders der zweiten 
Auflage, nur mit kleinen Änderungen und Ergänzun- 


gen versehen worden. Die Literaturhinweise wurden 


durch einige neuere Arbeiten ergänzt. 


Freiberg/Sa. _F. Rüns 


.J. Dieudonne, Foundations of Modern Analy- 
sis. (Pure and Applied Mathematics, Vol. X) XIV 
+ 361 S. New York 1960. Academic Press Inc. Preis 
geb. $ 8,50. : 


Das vorliegende Buch des berühmten französischen _ 
Mathematikers ist eine hervorragende Einführung in 
die heutige Analysis und Funktionalanalysis. Das in 
modernster Terminologie abgefaßte Werk setzt zwar 
keine besonderen Vorkenntnisse voraus, jedoch eine 
gewisse Fähigkeit zum abstrakten Denken. Eine 
große Anzahl von außerordentlich gut ausgewählten 
‚„‚Problemen‘‘, die hinter den meisten Paragraphen 
eingefügt sind, ergänzt den im Text dargebotenen 
Stoff und bietet einen Ausblick auf weitere wichtige 
und interessante Fragestellungen. Kap. I und Kap. II 
enthalten die wichtigsten Begriffe aus der Mengenlehre 
und aus der Theorie der reellen Zahlen. Kap. III, das 
sich eingehend mit den Eigenschaften der metrischen 
Räume beschäftigt, ist zugleich eine vorzügliche Ein- 
führung in ein Studium der analytischen Topologie, 
wie sie z. B. in den bekannten Büchern von KeELLEY 
und, BOURBARI dargestellt ist. Das nächste Kapitel, 
das sich nochmals mit dem Körper der reellen Zahlen 
befaßt, der aber nun im Anschluß an das vorherge- 
hende Kapitel als metrischer Raum betrachtet wird, 
schließt mit dem wichtigen TIETZE-URYsoHNschen 
Erweiterungssatz. Kapitel V, das auf normierte Räu- 
me und speziell auf BAnacH-Räume eingeht (Reihen 
in normierten Räumen, Teilräume und endliche Pro- 
dukte von normierten Räumen, äquivalente Normen, 
Räume von stetigen multilinearen Abbildungen, ab- 
geschlossene Hyperebenen und stetige Linearformen, 
Satz von F. Rınsz: Ein lokal kompakter normierter 
Raum ist stets von endlicher Dimension), wird jedem 
Leser, der sich besonders mit der allgemeinen Theorie 
der BanAcH-Räume und der topologischen Vektor- 
räume beschäftigen will, als Einführung von großem 
Nutzen sein. Auf einen Abschnitt über HILBERT-Räu- 
me folgen wichtige Ausführungen über Räume stetiger 
Funktionen (Approximationssatz von STONE-WEIER- 
STRASS, gleichstetige Mengen und der Satz von AscoLi1, 
Regelfunktionen). 50 Seiten sind dann dem natürlich 
völlig modern dargestellten Differential-Kalkül ge- 
widmet (der fundamentale Begriff der Differenzierbar- 
keit einer Abbildung aus einem BAnAcH-Raum in ei- 
nen anderen stützt sich auf die Approximation dieser 
Abbildung durch lineare Abbildungen), weitere 70 Sei- 
tenlegen ungemein elegant die Grundlagen der Theorie 
der analytischen Funktionen von mehreren Veränder- 
lichen dar. Dem Kap. X über Existenzsätze (Methode 
der sukzessiven Approximation, implizite Funktionen, 
Differentialgleichungen, Abhängiskeit der Lösung 
von einem Parameter usw.) schließt sich als Höhe- 
punkt des Buches das letzte Kapitel über die Grund- 
lagen der Spektraltheorie an (Spektrum eines stetigen 
Operators, kompakte Operatoren, Theorie von F. 
Rızsz usw.). Wie die Theorie auf bekannte Problem- 
stellungen der klassischen Analysis angewandt wer- 
den kann, wird besonders eingehend und instruktiv 
in den letzten beiden Paragraphen durch die Behand- 
lung der FreDHoLMschen Integralgleichung und des 
STURM-LiouviLLeschen Problems gezeigt. Ein ge- 
schickt ausgesuchtes Literaturverzeichnis beschließt 
das inhaltsreiche und äußerst empfehlenswerte Buch. 


Dresden M. LANDSBERG 
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B. Spain (Head of the Mathematies Department, 
Sir John Cass College), Analytical Quadries. IX + 
1358. m. 23 Abb. Oxford/London/New York/Paris 
1960. Pergamon Press. Preis geb. 30 s. net. 


Vorliegendes Werk ist vom Autor als eine kurz- 
gefaßte Einführung in die analytische Geometrie des 
dreidimensionalen euklidischen Raumes gedacht. Es 
wird dabei durchweg ein räumliches kartesisches 
Normalkoordinatensystem verwendet und auf die 
Hilfsmittel der Vektorrechnung verzichtet. Vom 
Leser werden nur die einfachsten Kenntnisse der 
analytischen Geometrie der Ebene, wie sie etwa in der 
Oberschule geboten werden, vorausgesetzt. Das 
1. Kapitel bringt als Einführung die kartesischen 
räumlichen Normalkoordinaten von Punkten, Ab- 
stands- und Teilungsformel, Richtungskosinus einer 
Geraden, Winkel zweier Geraden und die Projektion 
eines räumlichen Streckenzuges auf eine Gerade. Das 
zweite Kapitel ist der Ebene und der Geraden sowie 
den Grundaufgaben über Ebenen und Geraden ge- 
widmet. Die nächsten drei Kapitel behandeln Ku- 
geln, Kegel 2. Ordnung und Zylinder 2. Ordnung, 
Mittelpunktsflächen 2. Ordnung und das Paraboloid 
in einfachster Lage zum Koordinatensystem. Im 
6. Kapitel werden die Fernelemente durch Einfüh- 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


I. Szabö, Höhere technische Mechanik. 3. ver- 
besserte und erweiterte, Aufl. XII -+-5048. m. 
421 Abb. Berlin/Göttingen/Heidelberg 1960. Sprin- 
ger-Verlag. Preis geb. 31,50 DM. 


W. Traupel, Thermische Turbomaschinen. 
2. Band: Regelungsfragen, Festigkeit und dynami- 
sche Probleme. VIII + 4208. m. 459 Abb. Berlin/ 
Göttingen/Heidelberg 1960. Springer-Verlag. Preis 
geb. 61,50 DM. 


K. Klotter, Technische Schwingungslehre. 
Zweiter Band: Schwinger von mehreren Freiheits- 
graden (Mehrläufige Schwinger). 2. umgearbeitete 
und ergänzte Aufl. XVI-+483S. m. 296 Abb. 
Berlin/Göttingen/Heidelberg 1960. Springer Verlag. 
Preis geb. DM 58,50 DM. 
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% r x g ve. 
rung homogener Koordinaten erfaßt, der 2 5 
Kugelkreis in die Betrachtungen einbezogen und die 
Kreisschnitte der Flächen 2. Ordnung aufgezeigt. I 
Kapitel 8 bespricht die Eigenschaften einer allgemei- Fe 
nen Quadrik. Im 9. Kapitel werden die Invarianten 
einer quadratischen Form bei Schiebungen und 
Drehungen untersucht und seine Klassifikstion mit 
Reduktion von Quadriken durchgeführt (Hauptachsen- 
transformation). Im 11. Kapitel werden yomern 
schaften von Quadriken untersucht und konfo a 
Quadriken betrachtet. Sodann bringt Kapitel 2° 
lineare Systeme von Quadriken und deren schaf- 
ten. Im letzten (13.) Kapitel werden Ebenenkoordi- 
naten eingeführt und das Dualitätsgesetz im Raum 
studiert. Ein Anhang ist der Matrizenrechnung ge- 
widmet. 


Eine große Anzahl von Beispielen, deren Lösungen 
im 2. Anhang erscheinen, erleichtert dem Studieren- 
den die Erarbeitung des gebotenen Lehrstoffes. Für 
die Studenten der Technischen Hochschulen der 
ersten Semester ist daher diese auch didaktisch her- 
vorragende Einführung in die analytische Geometrie 
des Raumes sehr zu empfehlen. 
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Das Buch behandelt nach einer geschichtlichen Einleitung der Strom-Spannungs-Ma- 
trizen des Mehrtores (also die Leitwertmatrix, Widerstandsmatrix, Reikenparaltel, 
matrix, Parallelreihenmatrix und Kettenmatrix) seine Wellenmatrizen (scattering- 
matrix und na rAsenbunetree): die eine Veralleeneihernnn der klassischen Wellen- 


theorie des Vierpols darstellen. 


Die sehr einfache Strukturregel der Leitwertmatrix führt zu einem systematischen Be- 
rechnungsverfahren für die Vierpoleigenschaften beliebig komplizierter Schaltungen, 
das auch die Anwendung programmgesteuerter Rechenautomaten gestattet. Die 
weiterhin behandelte Transformationstheorie der Strom-Spannungs-Matrizen wird auf 
die Umwandlung der tee sie sowie auf die Fragen der Leitungsbildung in 
der Nebensprechtheorie angewendet. Als Anwendungsbeispiel der Vierpol-Wellenpara- 


metertheorie enthält das Buch schließlich die Theorie der Abzweigfilter und Brücken- 


filter, also der Grundfilter, der Zobelfilter und der Jaumannfilter. 
Bestellungen durch eine Buchhandlung erbeten 
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\ (Schriftenreihe der Institute für Me hei der De >. 
zu Berlin, Heft 9) En‘ 
1961. 149 Seiten — 5 Abbildungen — — gr. e— 200 


. Für die Theorie der BE Zahlen gibt BE 
‚weisen, die idealtheoretische, die auf Dedekind zu 
auf Grund der. p-adischen Methode von Hensel er 

- theorie von Ostrowski u, a. unterbaut wurde. 


In dieser Arbeit wird eine konsequente ;0 
(der Maximalordnungen und auch anderer Ordnungen) e 


e r r 
- . 


und durch eine Bewertungstheorie unterbatit, Be 


Das Buch ist so geschrieben, daß es. auch von Nichtkenmern der 
“ EReRiBe gelesen werden kann. 3 a 
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Theorie der Funktionen einer reellen Var 
Zweite ergänzte Auflage unter wissenschaftlicher Bearbeitung von K. Bag 


Übersetzung aus dem Russischen - 


(Mathematische Lehrbücher und Monographien, I. Abt., Band VI) 
1960. XII, 590 Seiten — 9 Abbildungen — gr. 8° — Ganzleinen DM 39,— 
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Elemente der Funktionalanalysis 
Zweite, unveränderte Auflage 


Übersetzung aus dem Russischen 


(Mathematische Lehrbücher und Monographien, I. Abt., Band VII) 
1960. XII, 256 Seiten — 4 Abbildungen — gr. 8° — Ganzleinen DM 25,— 
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